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Chapter 1 

Algebra, Trigonometry and 
Geometry 


A polynomial of degree n is an expression of the form 

P ( x ) = a n x n + a n -\x n ~ l + ... + a 2 x 2 + aix 1 + a 0 

n 

Sigma q ? ; = ai + a 2 + 03 + ••• + a n — 1 + a n 

i= 1 

Factorial n\ = n ( n — 1) (n — 2) ...3.2.1 
Finite series 

n 1 

^2 k = -n{n + 1) 
fc=i 

n 1 

A: 2 = -n(n + 1)(2 n + 1) 

' 6 
fc=i 

fe=i \fc=l / 

1.1 Indices and logarithms 

If x > 0 then we have the following properties: 

x° = l, x~ 1 = l/x, x m+n = x m x n , x m ~ n = x m /x n , x mn = (x m ) n x» = 
If a > 0 and a / 1 then 

log a 1 = 0, log a (.T _1 ) = log a (l/s) = - log a x 

log a x + log a y = log a (xy), m log a x = log a (x m ) 
loga loga y = logo (*/?/) 


1 



The logarithm and power functions are inverse functions, i.e. 
if x = log a y then y = a x 
and if x = a y then y = log a x 
Change of base log a x = log fe x/ log b a 

Logarithms to base e, i.e. log e x, are often written lnx. Such logarithms are 
called natural logarithms, e is the exponential constant given by 
e = l + V.+^.+i.+--- = 2-71828 .... 


1.2 Factors and roots of equations 


(a ± b) 2 = a 2 ± 2ab + b 2 
a 2 — b 2 = (a + b) ( a — b) 
a 3 ±b 3 = (a± b ) (a 2 =F ab + b 2 ) 

a n - b n = (a - b) (a"' 1 + a n ~ 2 b + a n ~ 3 b 2 + ... + ab n ~ 2 + b^ 1 ) 
Quadratic equation ax 2 + bx + c = 0 with roots a, (3 


a , (3 


-b±Vh 2 


4 ac 


/2a 


a + (3 = —b/a, a{3 = c/a 

Cubic equation ax 3 + bx 2 + cx + d = 0 with roots a, /3, 7 


a + (3 + 7 = — b/a, a(3 + f3y + 'ya = c/a, a(3'y = —d/a 


f (a) = 0 if and only if (x — a) is a factor of / (x). 


Remainder Theorem 

Suppose P(x) = a n x n + a n _ \x n ~ x + . . . + a\x + ao is a polynomial of degree 
n and that a is a root of the equation P(x) = 0. Then ( x — a) is a factor of 
P(x). 


1.3 Partial fractions 


If the degree of the polynomial f(x) is less than that of the denominator 
then: 


/ 0) 

{x + a) (x + b) ... 

/(g) 

( ax 2 + bx + c) (dx + e ) 

/(g) 

(ax + b) (cx + d) 2 


A B 

\ 7 + ••• 

x + a x + b 

Ax + B C 

ax 2 + bx + c dx + e 

ABC 

ax + b cx + d (cx T d) 2 
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1.4 Permutations and combinations 


The number of ways of selecting r objects from n objects with due regard 
to order is 

up = rc! 

r ( n — r)\ 

The number of ways of selecting r objects from n objects without regard to 
order is 

( n ) = — — (sometimes written n C r ) 

\ r J (n — r)\r\ 

Properties 

/ n \ _ / n \ /n+l\_/n\/ n A 

y r y y n — r J ’ y r J y r J y r — l J 

Binomial Theorem 

For any positive integer, n, 


(a + b) n 



a n-r b r 


1.5 Trigonometric functions 

Degrees and Radians 

360° = 2ir rad 1° = ^ rad 1 rad = (^)° w 57.296° 
cos nir = (— l) n , sinn7r = 0 

cos [(2n + 1) 7t/2] = 0, sin [(2 n + 1) 7r/2] = (—1)"' 

cos (tt/4) = sin (tt/4) = l/\/2 
cos (7r/3) = sin (7 t/ 6) = 1/2 
cos (7 t/ 6) = sin (tt/3) = \/3/2 

cos (tt/2 — 6) = sin 9 , sin (n/2 — 9) = cos 9 

Trigonometric identities 

cos -1 x = 7t/2 — sin -1 x, tan -1 x = tt/2 — tan -1 (1/x) 

sin 2 x + cos 2 x = 1, 

1 + tan 2 x = sec 2 x 
1 + cot 2 x = cosec 2 x 


sin (x + y) = sin x cos y + cos x sin y 
sin (x — y) = sin x cos y — cos x sin y 
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cos (x + y) = cos x cos y — sin x sin y 
cos (x — y) = cos x cos y + sin x sin y 


tan (x + y) = 
tan (x — y) = 


tan x + tan y 

1 — tan x tan y 
tan x — tan y 

1 + tan x tan y 


cos 2x = cos 2 x — sin 2 x = 2 cos 2 x — 1 = 1 — 2 sin 2 x 


sin 2x = 2 sin x cos x, tan 2x = 


2tanx 


1 — tan 2 x 

cos 3x = 4 cos 3 x — 3 cos x, sin 3x = 3 sin x — 4 sin 3 x 

x + y x — y 

sin x + sin y = 2 sin cos 

y 2 2 
x + y x — y 
sin x — sin y = 2 cos sm 


x + y 

cos x + cos y = 2 cos — r — cos 


2 ’ 
x-y 


2 2 
x + y x — y 
cos x — cos y = —2 sm sm 

& O O 


cos x cos y = \ [cos (x + y) + cos ( x — y)\ 
sin x sin y = \ [cos (x — y) — cos (x + y)\ 
sin x cos y = \ [sin (x + y) + sin (x — y)\ 

cos 2 x = \ (1 + cos 2x ) , sin 2 x = \ (1 — cos 2x) 

sin” 1 x ± sin” 1 y = sin -1 (xy/l — y 2 zb y\J 1 — x 2 ^ 

cos -1 x zb cos -1 y = cos^ 1 (xy =F \/l — x 2 -\/l — y 2 ^j 
tan -1 x zb tan -1 y = tan -1 [(x zb y) / (1 xy)\ 


x 


If t = tan — , then 
2 ’ 


2 1 


1 -t 2 


smx = 


cos x = 


tan x = 


2 1 


1 + t 2 ’ 1 + t 2 ’ 1-t 2 

cosx = (e? x + e~i x ) /2, sinx = ( e^ x — e~^ x ) /2 j 

De Moivre’s Theorem 

(cos x + j sin x) n = cos nx + j sin nx 
e? x = cos x + j sin x, e~ JX = cos x — j sin x 
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1.6 Hyperbolic functions 


coshx = ( e x + e x ) /2, sinh x = (e x — e x ) /2 
e x = coshx + sinh x, e~ x = coshx — sinhx 


tanh x = 


sinh x 


e — e 


0 2x 


- 1 


cosh x e x + e x 
sinh -1 x = In ^x + y/ x 2 + l^j 
cosh -1 x = ± In (x + V x 2 * — 1 j 


0 2x 


tanh 1 x = — In 
2 


1 + x 

1 — X 


1 - e -2 * 
1 + e~ 2x 


Relationship with trigonometric functions 


cosjx = coshx, coshjx = cosx 
sin jx = j sinhx, sinhjx = j sinx 


Osborne’s rule 

An identity involving hyperbolic functions may be obtained from the equiv- 
alent trigonometric identity by replacing the trigonometric functions with 
the corresponding hyperbolic functions and changing the sign of a product 
(or implied product) of two sines. 

e - g - 

cosh 2 x — sinh 2 x = 1, 1 — tanh 2 x = sech 2 x 


N.B. The Maclaurin series in powers of x for the hyperbolic functions may be 
found from the Maclaurin series for the corresponding trigonometric function 
by changing the sign of a product of two x’s 


e.g. 


2 4 

cosh x = l + — + — -) b 

2! 4! 


x 2n 

(2n)! 


+ ■■■ 


sinh 1 x = x — + — — — + . . . + (— 1) 


3x 5 5x 7 


(2n)! 


2n+l 


40 


112 


2 2n (n!) 2 (2n + 1) 


+ ... 


1.7 Geometry 
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Pythagoras’ theorem 


2 i j.2 2 

a + b = c 

Theorems 

The angle sum of a triangle is 180° or 7rrad. 

The sum of the interior angles of an n-sided polygon is (2 n — 4) x 90° or 
(n — 2)7 t rad. 

The exterior angle of a triangle is equal to the sum of the interior opposite 
angles. 

The angle subtended by a diameter at the circumference of a circle is 90°. 
The angle subtended by a chord at the centre of a circle is twice the angle 
subtended by the same chord in the opposite segment. 

The diagonals of a parallelogram bisect each other. 


1.8 Conic sections 

A conic section is the locus of a point that moves in a plane so that the ratio 
(the eccentricity) of its distance from a fixed point (the focus) in the plane, 
to its distance from a fixed line (the directrix), is a constant, e. 


1. Parabola (e = 1) ; focus at (a,0), directrix x = —a 
Cartesian equation: y 2 = 4ax 
Parametric equation: x = at 2 ,y = 2 at 


2. Ellipse (e < 1); foci at (±ae, 0), directrices x = ±a/e 

Major axis of length 2a, minor axis of length 2b 

y^ l 

Cartesian equation: — ^ = 1 with b = a (l — e 2 ) 2 

a z b z v ' 

Parametric equation: x = a cos 9, y = bsind 


3. Hyperbola (e > 1); foci at (=tae, 0), directrices at x = ±a/e 

y ^ 1 

Cartesian equation: — ^ — —r = 1 with b = a (e 2 — l) 2 
ci z b~ v ' 

Parametric equation: x = a sec 9, y = 6 tan# 

Rectangular hyperbola referred to its asymptotes as axes: xy = c 2 


4. Circle (e = 0) 

Cartesian equation x 2 + y 2 = a 2 
Parametric equation x = a cos 9,y = a sin 9 
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The polar equation for these three conic sections with the pole at a focus is 


- = 1 + e cos 9 
r 


The general conic 

The equation ax 2 + 2 hxy + by 2 + 2 gx + 2 fy + c = 0 represents: 
a circle if a = b and h = 0 

a pair of straight lines if af 2 + bg 2 + ch 2 = 2 fgli + abc 
a parabola if h 2 = ab 
an ellipse if h 2 < ab 
a hyperbola if h 2 > ab 
a rectangular hyperbola if a + b = 0 


1.9 Mensuration 

Circle , radius r: perimeter is 2nr. area is nr 2 

For a segment of angular width 9 (radians), arc length is rO and area is \r 2 9 
Ellipse , axes 2 a and 2b: perimeter is approximately 2it\J (a 2 + b 2 ) /2, 
area is irab 

Cylinder , radius r, height h: surface area is 27 t ?’ (h + r), volume is i rr 2 h 
Cone, base radius r, height h, slant height l : 

curved surface area is ir rl, volume is 7r?’ 2 /r/3 
Sphere, radius r: area is 47rr 2 , volume is 47rr 3 /3, 

area cut off by parallel planes distance h apart is 2irrh 
Triangle, ABC sides a,b,c: 
area of triangle is 

A = ^6csin7l = y/s (s — a) (s — b) (s — c), (Heron’s formula ) 
where 2 s = a + b + c 
Radius of circumcircle is R = abc/ 4 A 
Radius of inscribed circle is r = A/s 

a b c 

sine rule = — = — = 2 R 

sm A sin B sm C 

cosine rule a 2 = b 2 + c 2 — 2 be cos A 


sin (A/2) = yj {(s-b ) (s - c) /be} 
cos (A/2) = ^/{(s — a) /be} 

1.10 Complex numbers 

The algebra of complex numbers follows the same rules as that for real 
numbers with the addition of the unit imaginary number j which has the 
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property j 2 = — 1 or j = \/—l. 

Cartesian form z = x + jy where x is called the real part of z, Re(z), and 
y is called the imaginary part of z, Im(z). 

The Argand diagram is a geometric representation of the complex number 
z = x + jy, the point (x,y) represents z. 

Polar form z = r (cos 6 + j sin 0) where (r, 9) are the polar coordinates of 
(x,y) 

Exponential form z = re je 

Complex conjugate z = x — jy = r (cos 9 — j sin 9) = re~^ e 
Modulus \z\ = \J x 2 + y 2 , zz = r 2 = x 2 + y 2 
Argument (principal value) arg z = 9 , (— ir < 9 < ir). 

De Moivre ’s theorem (cos 9 + j sin 9) n = cos n9 + j sin n9 
Euler ’s formulae 

ei 6 = cos 9 + j sin 9, e~^ e = cos 9 — j sin 9 
cos 9 = \ (e je + e -J0 ) , sin 9 = jj (e J0 — e -je ) 

Complex roots If z = re^ e then the n complex roots of z are given by 

z V n = r i/n exp {j 6 + 2kn } k = 0, 1, 2, . . . , n - 1 
n 

These roots are equally spaced around the circle, radius r 1 /”, centred on the 
origin. 

Fundamental theorem of algebra 

A polynomial, P (z), of degree n given by 

P (z) = a n z n + a n -\z n ~ l H b a\z + a 0 ( a n / 0) 

can be factorized into n complex factors: 

P (z) = a n (z - zi) (z-z 2 )...(z- z n ) . 

The numbers z\, z 2 , • • • , z n are called the roots of the equation P (z) = 0. 

If the coefficients ao, a>i, . . . ,a n are all real then the complex roots occur in 
conjugate pairs. 



Complex variable 

If f(z ) = u(x,y) + jv(x,y) is an analytic function of the complex variable 
z = x + jy, then u and v satisfy the Cauchy- Riemann equations 

du dv du dv 

dx dy dy dx 

1.11 Inequalities 

Basic rules based on ‘greater than’ (the same rules apply to >,<,>,<) 

x>y<^x + a>y + a 

if p > 0 then x > y px > py 
if n < 0 then x > y nx < ny 
if a > b and x > y then a + x > b + y 
if x > y and y > z then x > z 

Vx + y <Vx + Vy 

if x,y > 0 then x < y x 2 < y 2 \fx < ^/y 
if x > 1 then x m > x n m > n 
if 0 < x < 1 then x m > x n <^> rn < n 
if x > 1 then x l ^ m > x 1 m > n 
Bernoulli’s inequality 

If x > — 1 then (1 + x) n > 1 + nx 
Arithmetic mean |( x -\- y ) Geometric mean y Gey 

\{x + y) > y/xy 

Triangle inequality |a; + 2/|<|a;| + |y| 

kl - \y\ < 11^1 - |y|| <\x-y\ 

Cauchy Schwarz inequality |u.v| < ||u||||v|| 

Minkowski inequality ||u + v|| < ||u|| + ||v|| 
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Chapter 2 


Determinants and Matrices 

2.1 Determinants 


2x2: , = ad -be 

c d 


3x3 : 

an 

021 

0.31 

a 1 2 

022 

032 

013 

023 

033 

022 

— On 

032 

023 

033 

021 

— 012 

031 

023 . 

+013 

033 

021 

031 

022 

032 



Oil 

012 

Oin 






n x n 

: A = 

021 

022 

02 n 








Onl 

On 2 

Oim 







The minor, otij , of the element is the (n — l)th order determinant formed 
from A by omitting the row and the column containing aij. 


The cofactor, A l} , of the element a ij is given by A tJ = (— lj b+3 a t j . 

The value of the n x n determinant is 

A = ctji Aj\ + a a A a + . . . + a m A m (expansion by ith row) 

= aij Aij + a-ijA-ij + . . . + a nj A nj (expansion by jth column) 
valid for i, j = 1 to n. 

Properties 

1 . The value of the determinant is unchanged if its rows and columns are 
interchanged. 

2. The value of a determinant is unchanged if a multiple of any row (or 
column) is added to any other row (or column). 

3. If the elements of two rows (or columns) are proportional, i.e. linearly 
dependent, then the value of the determinant is zero. 
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4. If the elements of two rows (or columns) are interchanged then the 
resulting determinant has the same numerical value but its sign is 
changed. 

5. If the elements of a row (or column) are all multiplied by a constant 
k, then the value of the determinant is also multiplied by k. 

Cramer’s rule (efficient for only two equations in two unknowns or three 
equations in three unknowns). 


ail.Tl + 012^2 + 013X3 

= hi 


Oil 

012 

013 

021^1 + <122X2 + 023^3 

= h 2 

A = 

021 

022 

023 

031^1 + 032X2 + 033X3 

= hi 


031 

032 

033 



hi 

012 

013 


Oil 

hi 

013 


On 

012 

hi 

Ai = 

h 2 

022 

023 

! ^2 = 

021 

h 2 

023 

! ^3 = 

021 

022 

h 2 


h 3 

032 

033 


031 

h 3 

033 


031 

032 

h 3 


The solution is 


xi = Ai/A, x 2 = A 2 /A, x 3 = A 3 /A 


2.2 Matrices 

The m x n matrix is written as 


A = 


or A = \a . 


■ij\ 


on 

012 • • 

din 

021 

022 • • 

02 n 

Oml 

Om2 

dmn 

is n 

x m 



A is called a symmetric matrix if A 7 = A 

If A and B are of the same order then qA + /IB = [aa,j + /3bij] 
If A is m x r and B is r x n then 


AB = 


^ ^ Uikbk : j 


Lfc=i 


is m x n 


Unit matrix or identity matrix I = [5ij\ where the Kronecker 5 is given by 


$ij — 


0 , i^j 

1 , * = 3 
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Null matrix or zero matrix 0 = [0] 

For a square matrix A if det A ^ 0 then the inverse matrix , A -1 , exists 
and is given by A -1 = dc * A adj A, where adj A is the transposed matrix of 
cofactors. 

The inverse matrix has the property AA _1 = A~ X A = I 
If det A = 0 then A is said to be singular. 

A is called an orthogonal matrix if A 1 2 = A -1 
For a 2 x 2 matrix 


a b 

~ x _ 1 

d 

-b ' 

c d 

ad — be 

—c 

a 


For an n x n matrix 

det A = det A r and det ( kA ) = k n det A 

The rank of a matrix, A, is the largest non-zero sub-determinant of A. 

A set of vectors {vi, V 2 , . . . , v„,} for a finite dimensional space is said to be 
linearly dependent if there exist scalars aq, a 2 , . . . , a n , not all zero such that 

aqvi + a 2 v 2 + . . . + a n v n = 0 

otherwise the set is linearly independent. 

Suppose that {vi, v 2 , . . . , v n } is an arbitrary basis for a finite dimensional 
inner product space. Then the orthogonal basis {wj, w 2 , . . . , w n } may be 
obtained by the Gram Schmidt process 

< V n , W[ > < v n ,w 2 > < v n ,w n _i > 

w n = V n n jrq Wi n j-q W 2 , . . . , -5 W n _i 

ll w ilr ll w 2 || 1 1 w n — 1 |p 

2.3 Systems of equations 

For a system of m linear equations in n unknowns we write in matrix form 


Ax = b 


The explicit form is 


anxi + ai 2 x 2 + 

~\~CL\ n X n 

= for 

a 21 X 1 + a 22 x 2 + 

-\-&2n%n 

= b 2 

®ml®l T O m2 X 2 - 1- 

~\-Q"mn%n 

= b m 


1. If rank [A,b] > rankA then there is no solution. 

2. If rank [A,b] = rankA = n then there is a unique solution. 
In the case m = n, the solution is given by x = A *b. 
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3. If rank [A,b] = rank A < n then there is an infinite number of solu- 
tions. 

It follows from 2. and 3. that, for a square matrix A, the homogeneous set 
of equations Ax = 0 has a non-trivial solution if and only if detA = 0. 

2.4 Eigenvalues and eigenvectors 

For a square matrix A, the eigenvalues A ; and corresponding eigenvectors 
Xj are given by 

Ax, = A jXj 

Thus the eigenvalues satisfy det [A — AI] = 0 

Properties 

Matrix pA A 1 A -1 A k A + ql [A+pI] -1 

Eigenvalues p\, A * 1/A, X k A i + q 1/ (A i + p) 

In each case the corresponding eigenvectors are x,. 

The trace of the matrix A, is given by 

trA = an + o 2 2 + ■ • • + ^nn 

A i = trA, n^ = ^ 

Suppose that A is an eigenvalue of A then A = 0 iff det A = 0. 
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Chapter 3 

Vector Algebra 


3.1 Vector addition 



Triangle law 
Component form 


a + b 
a + b 


c 

(ai + 61)1 + (a 2 + b 2 ) j + (a 3 + b 3 ) k 


3.2 Vector products 



0 is the angle between the two vectors a and b, n is the unit vector perpen- 
dicular to both a and b in the sense that a, b and n, in that order, form a 
right-handed system. 

a.b = ab cos 0 = b.a, axb = a&sin0n 
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In component form 


a = aii + a 2 j + a 3 k, b = 6 x i + b 2 j + 63k 


a.b = aifti + a 2 b 2 + a 3 b 3 


axb = 


i J k 

a x a 2 a 3 
b\ b 2 b 3 


= — bxa 


Magnitude: |a| = a = yaf + a| + a| 

Unit vector in the direction of the vector a is a = - — -a = 

a 


a. (bxc) 


(axb) .c = 


ai cl 2 a 3 

61 b 2 b 3 

ci c 2 c 3 


[a b c] 


a x (bxc) = (a.c) b — (a.b) c 
(axb) x c = (a.c) b — (b.c) a 


(axb) . (cxd) 


a. c a.d 

b. c b.d 


(axb) x (cxd) = [a b c] c — [a b c] d 


3.3 Polar coordinates in two and three dimensions 

Plane polar coordinates (r, 6) 

x = r cos 0 , y = r sin 6 
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Cylindrical polar coordinates (R, (j), z) 

x = R cos 4>, y = R sin (j) 



Relationships with Cartesian unit vectors: 

R = cos 4>i + sin , 0 = — sin cj)i + cos </>j 

Spherical polar coordinates (r, 0 , (f>) 

x = r sin 0 cos 4>, y = r sin 9 sin </>, z = r cos 9 



Relationship with cylindrical unit vectors: 

r = sin 0R + cos 0k, 6 = — cos 0R + sin 0k 
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Chapter 4 

Calculus 


4.1 Rules for manipulation of derivatives and integrals 


Differentiation 

r . . , d , . du , dv 

Linearity ride — ( au + bv) = a— — \- b— 
dx dx dx 


Product rule 


d , . du dv 

— (uv) = —v + u— 
dx dx dx 


Quotient rule 
Leibniz rule 


d /u\ f du dv\ 


dx \v 


dx 


dx J 


h 7 


D n (uv) = uD n v + ( i ) ( Du ) D n ~ 1 v + ( 9 ) ( Du ) ^‘" 2 « + • • • + (■ D n u ) v 

( D r u ) D n ~ r v, where D = — 

dx 

Chain rule for ordinary differentiation (function of a function rule) 

Leibniz notation 

y = y (x) and x = x (it) ; 
dy dy dx 

du dx du 

Function notation 

[F (g (x))]' = F' (g (x)) g' (x) 

Chain rule for partial differentiation 

F (it, v) = f (x, y) with x = x (u,v) ,y = y (u, v) 

dF df dx Of dy dF df dx df dy 

du dx du dy du ’ dv dx dv dy dv 
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Total differential 


df 


df df 

-f-dx + -f-dy 
ox oy 


Integration 

Integration by parts 


fdv [du 

/ u—dx = uv — / v—dx 
J dx J dx 


Integration by substitution 


f (x) dx 



dx 

u ) —du where F (u) 
du 


f (x («)) 


f r (x) 

J f O) 

/'(*) [/(s)] n d® = 


dx = In |/(x)| 


[/ (*)] 


71+ 1 


n + 1 


n / — 1 


Differentiation of integrals 

rv(x) 

If / (x) = / F(x, then 


J u(pc) 


dl_ 

dx 


/ \ \ dv Tn , . ..du 

I (x, vlx))— I (x, u(x ))- — h 

dx dx 


r v{x) Qjp 


/ u( x) 


dx 


(x, t)dt 
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4.2 Standard derivatives and integrals 

Table of derivatives and integrals 1 


df jdx f{x) F(x ) = J f(x)dx 

(add arbitrary constant where necessary) 


nx 


71—1 


1 1 

x 2 x 


a x In a a x 


™n+ 1 

In \x\ 


In a 


1 

X 

lnx 

x(lnx — 1) 

1 

log b x 

x(log b x - log fe e) 

x In b 

cos X 

sinx 

— cosx 

— since 

cosx 

sinx 

sec 2 x 

tanx 

In secx 

sec x tan x 

secx 

J In secx + tan x| 

\ In tan(x/2 + 7t/4)| 

—cosecx cot x 

cosecx 

J In cosecx — cot x| 

\ In tan(x/2)| 

— cosec 2 x 

cot X 

In | sinx| 

1 

sin -1 x 

x sin” 1 x -(- v 1 — x 2 

\J\ — X 2 

-1 

cos” 1 X 

x cos” 1 x — V 1 — x 2 

a/1 — x 2 

1 

tan” 1 x 

xtan” 1 x — - ln(l + x 

1 + X 2 

coshx 

sinh x 

coshx 

sinh x 

cosh x 

sinh x 

sech 2 x 

tanh x 

J In cosh x 
| ln(e x + e _a; ) 


(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 

(4.17) 

(4.18) 


continued. . . 
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Table of derivatives and integrals 2 


df/dx /(x) 


F(x ) = J f(x)dx 

(add arbitrary constant where necessary) 


— sechx t anh x sechx 

—cosechx coth x cosechx 

— cosech 2 x coth x 

sinh -1 x 
cosh -1 x 
tanh -1 x 

1 


V 1 + x 2 

1 

\Jx 2 — 1 

1 

1 — X 2 


X 


(a 2 — x 2 ) 

3 

2 

\J a 2 — x 2 

— 2x 


1 

(x 2 + a 2 

j 2 

x 2 + a 2 

a 2 — 2x 2 


1 

x 2 (x 2 — a 2 ) 

T 

2 

x\/.t 2 — a 2 

—x 


1 

(a 2 + x 2 ) 

T 

2 

Va 2 + ic 2 

—x 


1 

(x 2 — a 2 ) 

T 

2 

a/.t 2 — a 2 

2x 


1 

(a 2 — x 2 


(a 2 — x 2 ) 


2 tan - 1 (tanh x /2) 

2 tan -1 
tan -1 (sinhx) 

In | tanh(x/2)| 

— In | cosechx + coth x| 
In \(e x — l)/(e x + 1)| 

In | sinhxj 
In \e x — e~ x \ 


x sinh 1 x — yl + x^ 


x cosh 1 x — \J x 2 — 1 


x tanh 1 x — - ln(l — x 2 ) 


sin 


-l 


— cos 


-1 x 
■ tan — 
a a 


«x 


-l * 
- sec — 
a a 


-l 


x 


sinh 

a 

In |x + \J x 2 + a 2 1 

cosh -1 — 
a 

In |x + \/ x 2 — a 2 1 

1 , _i x 

— tanh — 
a a 


— in 
L 2a 


a + x 


a — x 


(4.19) 

(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 

(4.28) 

(4.29) 

(4.30) 


continued. . . 
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Table of derivatives and integrals 3 


df/dx f(x) 


X 

Vx 2 ± a 2 

Vx 2 ± a 2 

—x 

\J a 2 — x 2 

Va 2 — x 2 

e ax ( a cos bx — b sin bx) 

e ax cos bx 

e ax ( a sin bx + b cos bx) 

e ax sin bx 


F(x) = J f(x)dx 

(add arbitrary constant where necessary) 

2 

^ \/ x 2 ± a 2 ± In | x + \J x 2 ± a 2 | (4.31) 

x /—z a 2 x . , 

— \J a 2 — x 2 + — sin — (4.32) 

Z z cl 

e ax 

-^(acosbx + bs'mbx) (4.33) 

a 2 + b 2 K ' 

e ax 

-r, -^(asinftx — 6 cos bx) (4.34) 

a 2 + b 2 y ' 


4.3 Definite integrals 

Wallis’s formulae (reduction formulae also hold if upper limit is 7 r or 27 t. 


S n = / 2 sin" 6d0 = - — — 5'n — 2 

Jo* n _ 

C n = [ 2 cos n Odd = - — -C n _ 2 . 

Jo n 


t f 2 ■ mn n n in Ul ~ 1 r U-l 

Irn r, = I Sill 0 COS udu = i m -2,n = 4m.r7.-2. 


m + n 


m, + n 


'■m,n— 2* 


_ (n- l)(n-3)(n-5)... 
n " n(n — 2)(n — 4) . . . ’ 

r _ (m — 1 )(m — 3) . . . (n — l)(n — 3) . . 

4 m, n = yry— 7—777— —— o\ / , 7 \ F 


(m + n)(m + n — 2 )(m + n — 4) 


. 1 ^ n even 

where p = < 7 

1 1 n odd 


_ 1 | both m and n even 
^ '1 otherwise 


For all integer values of m and n: 


sin rnO cos nddO = 0 


sin nr# sin nddO = 


cos m6 cos nddO = ir6 n 


(8 rnn is the Kronecker 5 see page 11) 


m, n / 0 


21 




/ oo 

e~ x2 dx = v/tt; 

-OO 

dx = 7re _ l a l; 


/ 


cos ax 
72 1 


— OO 1 “I - ^ 

Error function 


f e-^ ax2+bx) dx = ^exp(6 2 /4a) a > 0 

/ 


sm ax 
— oo 4 “l - x 


dx = 0 


2 f x 2 

erf (re) = — ^ / e _< dt 

Jo 

Complementary error function 

2 

erfc(x') = 1 — erf(x) = — p= / 

V7T A 


e <2 (Zt 


erf(O) = erfc(oo) = 0, erf(oo) = erfc(0) = 1 
For the Chebyshev polynomials T n (x), T m (x) (see page 57) 


f 1 T n (x)T. m (x) 
/- 1 \/l — x 2 


0 n / m 

dx = { 7 r n = m = 0 

7r/2 n = m / 0 


Gamma function 

/»oo 

T(t) = / e~ x x t ~ 1 cix t / 0, — 1, — 2, . . . 

r(t) = ( °- i)r(t - 1 ) 

= (t - I)(t - 2)(t - 3) . . . r(i - [t]) 

r(±) = v^ _ 

r(n + 1) = n! when n is an integer 
5eta function 

B(s,t) = f x s_1 (i — x) t_1 dx = r(s)r(t)/r(s + 1) 
Jo 


Legendre polynomial 


P “ w = h i lx + (* 2 - !)” 

P 0 (®) = 1, Pi(x) = X, P 2 {x) = (3x 2 - l)/2, P 3 (x) = (5x 3 - 3x)/2, . . . 

Generating function 


r* 27T 


1 d" 


(1 — 2xh + h 2 ) 2 = ^ P n {x)h r 


n = 0 
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P n (x) is a solution of Legendre ’s equation (1 — x 2 )y" — 2 xy' + n(n + 1 )y = 0. 
A second, linearly independent, solution is the Legendre function of the sec- 
ond kind Q n (x ) which is not a polynomial. 

Bessel function 

/* "2" 

Jrt(x) = ; i— / cosfxsin#) cos 2 ” 9d6. 

2- 1 v^Fr(n + i) 7o 

J n (x) is a solution of Bessel’s equation x 2 y" + xy' + (x 2 — n 2 )y = 0 
A second, linearly independent, solution, Y n (x), is the Bessel function of the 
second kind. 

The modified Bessel function , I n (x), satisfies the equation 

x 2 y" + xy' — (x 2 + n 2 )y = 0 and a second linearly independent solution is 

K n (x), the modified Bessel function of the second kind. 


4.4 Radius of curvature of a curve 


Intrinsic coordinates: p = 
Cartesian coordinates: p = 
Parametric coordinates: p 


— . Curvature: k 

dip 

[1 + {y'?f /2 
\y"\ 

_ ((x) 2 + (y) 2 ) 3/2 
I xy ~ xy | 


1 

P 


4.5 Stationary points 

Functions of one variable 

The function /(x) has stationary points given by /'(x) = 0. The stationary 
point, (xo,/(xo)), is a local maximum if f"(x 0 ) < 0 or a local minimum if 
f"(x 0 ) >0. If f"(x 0 ) = 0 then the second derivative test is inconclusive 
and we consider the sign of /'(x ± e) where e is small and positive. If 
f"{x 0 ) = 0 and f{x 0 ± e) are both positive or both negative then xo is a 
point of inflection. 

Functions of more than one variable 

Stationary points occur when the first partial derivatives vanish. For a 
function of two variables, /(x, y), stationary points are given by f x = f y = 0. 
The Hessian is given by 


H(x,y) 


fxx 

fyx 


fxy 

fyy 


The stationary point is: 
a maximum if f xx < 0, f yy < 0 and H > 0, 
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a minimum if f xx > 0, f yy > 0 and H > 0, 
a saddle point if H < 0. 


Lagrange multipliers 

To find the stationary point of f(x,y) subject to the constraint g(x,y) = 0, 
find the stationary points of the function 

</>0, y; A) = f{x, y) + A g(x, y). 


4.6 Limits and series 

sinx „ tan x _ / x\ n _ 

lim = 1, Inn = 1, Inn 1 H = e . 

x — >0 X a;^0 X n—KX \ nJ 

lim x Q lnx = 0, lim x~ a lnx = 0, lim x a e~ x = 0, (in all cases a > 0). 

>0 x — >oo x—>oo 

De L’Hopital’s rule 

If f(a) = g(a ) = 0 , then lim ^7—^ = lim ^77— r, provided g'(a) / 0 . 


y(*) 


g'(x) 


Taylor’s theorem 

f{x) = f(a) + {x - a)f'(a) + f"(a ) + . . . + — — ^-/ (n) («) + 


2 ! 


n! 


(x-a) n+1 f+1 x 

where e n = — P '(c) for some c £ a, x). 

(n + 1)! w v ; 

Taylor’s series for a function of one variable 


/(x) = /(a) + (x - a)f'(a) + 
or, writing x = a + h, 


(x — a)^ 

2 ! 


f"( a ) + • • • + 


(x — a) 1 


n\ 




h 2 Un 

/(a + h ) = /(a) + hf\a) + + • • • + ^j/^( a ) + • • • 

Maclaurin series (put a = 0, h = x in Taylor’s series) 

2 

/(x) = /( 0) + x/'(0) + | r r ( 0) + . . . + |y/ (n) (0) + • • • 

Power series 

OO 

If the power series ao + aix + a 2 X 2 + . . . + a n x n + . . . = a n x n is convergent 

n = 0 

for |x| < R then R is called the radius of convergence of the series. 

R = lim 

«— > 00 a„ 
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-1 < X < 1 


Some useful series: 

/y>^ /y»0 ry* ^ /•»»*' 

ln(l+x) = X- — + y - — + + (-l) n+1 — 

/y»2 /vi3 /y»4 /yiH 

_ / V if* il/ »X/ 

ln(l-x) = -x-- •••, -1<x<1 

2 3 4 n 

. n(n — 1) o , n(n — l)(n — 2) o 

(1 ± x) n = 1 ± nx + — - — j — - x 2 ± — y -x 3 + ... - 1 < x < 1 

[This is the binomial series. If n is a positive integer the series terminates 
after n + 1 terms, otherwise it converges if and only if |x| < 1.] 


(IT®) 


-i _ 


= 1 ± x + x z ±x 6 + . . . + x 2n ± X 


^271+1 


+ .... 


Ixl < 1 


e = 


sm x = 


cos x = 


tanx = 


sin 1 x 


tan 1 x 


x 2 x 3 x 4 x n 

i + x + — + — + — + all x 

x 3 x 5 x 7 (— l) n x 2n+1 

' T “¥ + ¥ _ 7! + " + (2n+l)! + '"’ 

x 2 x 4 x 6 (— l) n x 2n 

1 1 1-... + - — b 

2! 4! 6! (2n)! 

x 3 2x 5 x 7 


all x 
all x 


7T 


X + T + T^ + 315 + ---’ N< 2 

n -l 

COS X 

2 

x 3 3x 5 5x 7 (2 n)! x 2n+1 

• T+ y + l0' + Tl2 + ”- + 22n(n!) 2 (2n + l) + '”’ N<1 

...3 ^5 ^2n+l 


7 y» ^ /v> v /y»^“ 

I -T + T + '" + (- 1 >” +1 7f-m 

3 5 (2n + 1) 


x < 1 


To obtain the series for the corresponding hyperbolic functions, see the note 
following Osborne’s rule on page 5. 


Taylor’s series for a function of two variables 

f(a + h,b + k) = f(a,b) + {hf x (a + b) + kf y (a,b)} 


n = 0 


dx dy 

4.7 Multiple integration 


+^{h 2 fxx(a, b ) + 2 hkf xy {a, b ) + k 2 f yy (a , 6 )} + . . . 

£ ¥ ('*/¥ + *!r) 


f(x,y)dxdy= / f(x(u,v),y(u,v)) 

R J JR' 


J 


x,y 

U , X 


dudv 


where the region in the xy-plane is mapped to the region R' in the mi- 
plane and the Jacobian of the transformation is given by 


J 


( - 

Xu 

Uu 

\u,vj 

Xy 

Uv 
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Plane polar coordinates 


J = r, so that dA = dxdy = rdrdO 



Volume and surface integrals 

Cylindrical polar coordinates 

dS = Rd(j)dz and dV = RdRdcf)dz 
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Spherical polar coordinates 


dS = r 2 sin Odddcj) and dV = r 2 sin OdrdOdq i> 



Line integrals 

Green's theorem in the plane 

L iPdx+ Q iv) =J L(^-%) dxdti 

where C is the closed curve bounding the region R. 

/ (Pdx + Qdy ) is independent of the path if and only if — — 

Jc dx 

The line integral of F is independent of the path C joining A to . 
if there exists a potential function (/> such that F = grade/*, or 
if and only if cur IF = 0. In this case 



f' 


d4> = 4 >b — 4>a 


4.8 Applications of integration 

Arc length of a curve 

r b 

Cartesian form: / \/l + (i/) 2 dx 
J a 

f° 2 !— 

Polar form: / -\J r 2 + ( dr/d9) 2 d6 


> 0 1 


rt2 


Parametric form: / \/ x 1 + y 2 dt 


_ dP 
dy ' 

3, if and only 
equivalently, 
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Solid of revolution 

The surface area of the solid of revolution obtained by one complete revo- 
lution of the curve y = f(x) about the x-axis is given by 

r b 

27 r / yyj 1 + ( y') 2 dx 
J a 

The volume of the solid of revolution obtained by one complete revolution 
of the curve y = /(x) about the x-axis is given by 

f b 

7 r / y 2 dx 
J a 

Centroid 

The centroid of the volume of the solid of revolution above is at (x, 0 ) where 



The centroid of the area between the curve y = /(x), the x— axis and the 
lines x = a, x = b is given by 


x = 



V = 





Mean value 

The mean value of /(x) between x = a, x = b is 


1 f b 

7 / f{x)dx 

b~ a J a 

Root mean square value 

The RMS of the function /(x) between x = a,x = b is given by 


1 

b — a 




1 

2 


Moments of area and inertia 

The first and second moments of a plane area A about an axis are given 
respectively by 


rdA and 


r 2 dA 


where r is the distance from the axis of the element dA. 
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The moment of inertia, /, of a body, of density p and volume V, about 
an axis is given by 

[ r 2 pdV 

Jv 

where r is the distance from the axis of the element dV. 


Parallel axes theorem 

If Iq is the moment of inertia about an axis through the centroid and I is 
the moment of inertia about a parallel axis distance d away, then 

I = Iq + rn.d 2 

Table of moments of inertia 1 


Uniform Body 

Axis 

M.of I. 

mass m 



1. Bar length 2a 

Perpendicular to bar 
through one end 

ma 2 

3 


Perpendicular to bar 

4 rrm 2 


through centroid 

3 

2. Rectangular lamina 

Parallel to side 26 

ma 2 

sides 2a and 26 

through centroid 

3 


Perpendicular to plane 

m(a 2 + 6 2 ) 


through centroid 

3 

3. Rectangular solid edges 

Perpendicular to face 4a6 

m(a 2 + 6 2 ) 

2a and 26, depth 2d 

through centroid 

3 


Coinciding with edge 2d 

4m (a 2 + 6 2 ) 

3 



2m?’ 2 

4. Solid sphere radius r 

Diameter 

5 



7 mr 2 


Tangent 

5 

5. Spherical shell radii 

Diameter 

2m. (R 5 — r 5 ) 

R and r 

3 {R 3 - r 3 ) 

m{7R 5 — 5 R 2 r 3 — 2 r 5 ) 

Tangent 


5(R 3 — r 3 ) 


continued. . . 
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Table of moments of inertia 2 


Uniform Body 

Axis 

M.of I. 

mass m 



6. Disc radius r 


mr 2 

Diameter 



4 


Perpendicular to disc 

mr z 


through centroid 

2 

7. Annular ring radii 

Diameter 

m(R 2 + r 2 ) 

R and r 

4 


Perpendicular to plane of 

m(R 2 + r 2 ) 


ring through centroid 

2 

8. Solid cylinder radius 

Axis 

mr 2 

r length 2a 

2 


Perpendicular to axis 

fa 2 r 2 \ 

“(t + tJ 


through centroid 

9. Cylindrical shell 
radii R 

Axis 

m(R 2 + r 2 ) 

o 

and r length 2a 


Z 


Perpendicular to axis 

( a 2 R 2 + r 2 \ 


through centroid 

”1.3 + 4 ) 


For entries 2, 6 and 7 the second moment of area J is given by J = IA/m 
where A is the area. 
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4.9 Kinematics and dynamics 

Kinematics 


Position vector 


r 


xi + yj + zh (Cartesian coordinates) 
r? (Plane polar coordinates) 


Velocity vector 


r 


xi + yj + ik 
rr + 

•S't 


(Cartesian coordinates) 
(Plane polar coordinates) 
(Intrinsic coordinates) 


Acceleration vector 


r 


xi + yj + z k 
(r — r0 2 )r + 

s 2 ^ 

st H n 

P 


ld_ 

r dt 


(r 2 d)0 


(Cartesian coordinates) 
(Plane polar coordinates) 

(Intrinsic coordinates) 


For uniform motion with angular velocity w in a circle of radius a, speed 
is v = au and acceleration is v 2 /a = aw 2 directed towards the centre. 


Dynamics 

Newton’s second law for a particle of mass m moving under the influence of 
an external force F 

F = mr 

For a rigid body moving about an axis with angular velocity lu and external 
torque P we have 

T = Id) 

where I is the moment of inertia about the axis. 
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Chapter 5 


Ordinary Differential 
Equations 


5.1 First order equations 


^ = F{x,y) or y = F{x,y) 

Variables separable 

^ = f{x)g{y) (. g(y ) j L - 0) 

General Solution is 



f{x)dx + constant 


Linear equation (integrating factor method) 

^ + P(x)y = Q(x) 

Integrating factor: g(x) = exp (f P(x)dx) 

Then —{g(x)y) = fi(x)Q(x ) and the 
dx 


General Solution is 

n(x)y = J /i(x)Q(x)dx + constant 

Initial-value problem 

Find the general solution of 

= F(x, y) subject to y(x Q ) = y 0 

i.e. solve the differential equation containing one arbitrary constant and 
determine the value of the constant to satisfy the initial condition that 
y = t/o when x = xq. 
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5.2 Second order equations 

The general linear equation with constant coefficients 

d 2 y . ,dy //,,/, ,, x 

a da? + b dx + Cy = ^ ° r ay +by + Cy = 

Complementary function y c 

Solution of the associated homogeneous equation 

ay" + by' + cy = 0 (a, b, c real) 

Auxiliary equation a\ 2 + bX + c = 0 
Three cases: 

Two real distinct roots, Ai, A 2 

Linearly independent solutions e Xlx and e X2X \ 

Vc = C ie x ' x + C 2 e^ x 


One real repeated root, Ai 

Linearly independent solutions e XlX and xe XlX \ 

y c = e Xlx (C 1 +C 2 x) 

Complex conjugate roots, a ± j/3 

Linearly independent solutions e ax cos fix and e ax sin/3x; 

y c = e ax (Ci cos (3x + C 2 sm/3x) 

Particular integral y. p 

Any solution of the non-homogeneous equation 

ay" + by' + cy = f(x ) 

Method of undetermined coefficients 

Choose a suitable trial function for y p according to the right-hand-side, f(x). 


/(*) 

trial function 

constant 

constant 

px + q 

lx + rn 

px 2 + qx + r 

lx 2 + mx + n 

polynomial degree n 

a n x n + a n _ix" -1 + . . . + aix 1 + a 0 

e “ 


p COS LUX 

m cos lux + n sin lux 

q sin lux 

m cos lux + n sin lux 

Re[(p — jq)e^ x ] 

R e[ze ju}X ] 


Beware: If the rhs includes one of the complementary function linearly 
independent solutions then choose x times the appropriate trial function. 


33 



General Solution The sum of the complementary function and a particular 
integral y = y c + y p . 

Initial-value problem 

ay" + by' + cy = f(x) subject to y{x) = y 0l y'(x 0 ) = y' 0 

Find the general solution containing two arbitrary constants and then use 
the initial conditions to determine the values of the constants. 


5.3 Higher order equations 

The general linear equation with constant coefhcents 

anV^ + a n _iy (n-1) + . . . + a 2 y" + aiy + a 0 y = f(x) 

The complementary function, y c , is found in a similar manner to that on 
page 33 using the auxiliary equation 


a n \ n + a n _iA fi 1 + . . . + aiA + fflo — 0 

A particular integral, y p , is also found in a similar manner to that on 
page 33 using suitable trial functions. 

The general solution is y = y c + y p . 
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Chapter 6 

Fourier Series 


Full-range series 

/(x) is defined for —l<x<l. 

The Fourier series for /(x) is given by: 

oo oo 

1 U7TX \ ^ 7 • U7rx 

F(x) = -a 0 + 2_^a n cos — h 2^ K sin — — 


where 


a ° ~ 1 


f(x)dx, 


’-l 


rl 


I ff \ nvx 
a n = j j /(xjcos —j—dx-, 


bn. — 


f l f( x • rmx 
j j[x J sm — - — ax 


F(x) extends /(x) periodically with period 21. 

F[x ) = /(x) at all points, —l<x<l, at which / is continuous. If / is 
discontinuous at xq, then F(x o) = {f(x o+) + f(x q— )}/2. 


Half-range series 

/(x) is defined for 0 < x < l. 

The Fourier half-range cosine series is: 

"I CXD nl 

1 U7TX Z / T17TX 

F c {x) = T^ap + )> cos — — ; a n = j /(x)cos — — dx. 


F c (x) extends /(x) periodically as an even function with period 21. 
The Fourier half-range sine series is 


OO 

F s {x) = ^ b n sin 

71 — 1 


n7rx 

“T 


brt. — . 


, n . nnx 
y(x) sm — - — ax. 


F s (x) extends /(x) periodically as an odd function with period 21. 

The Fourier series of some frequently occurring functions are given on page 47. 
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Complex Fourier series f(x) is defined for —l<x<l. 

The complex Fourier series for f(x) is given by 

OO 

F{x)= ]T c n e^ xjl 

n =— oo 

where 

i nl ( (&n Tl > 0 

Cn = 9 / / f( x ) e ~ Jnnx/ldx = { «o/2 n = 0 

l (a- n +jb- n )/2 n< 0 

The constant term, i.e. |ao or co, is the average value of f(x) over the 
interval (—1,1). 


Parseval’s theorem 

2 


1 °° 1 

[f(x)] 2 dx = -0§ + Y n(«n + bl ) 
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Chapter 7 

Vector Calculus 

7.1 grad, div and curl 


d ? d - d 
Nabla operator V = i— — b j — — b k— 
ox dy dz 

grad0 = V0 = —i + — j + — k 
ox Oy Oz 

dlvF = V.F = ^i + "k + BF > 


dx dy 


dz 


curlF = V x F = 
,dF 3 dF 2 


dy 


dz 


i + 


i j k 

d/dx d/dy d/dz 
F\ F-2 F : > 

dF\ dF :i ~ BF 2 dF ] 


dz dx 


■)j + (■ 


dx dy 


)k 


, , < 90 . 

grad cp = — n 
dn 


where n is the unit vector normal to the surface 0 = constant. 


Directional derivative. The rate of change of a scalar field, 0, in the direction 
of the unit vector n is given by 

<9 0 , , . 

— = grad (p.n 
dn 

If 6V is volume bounded by the surface 5S then 

1 f 

grad 0 = lim — 7 / ddS 
sv^o dF Jss r 

div A = lim —— / A.dS 

sv^o 5V Jss 
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curlA = lim 


1 


< 5 V ^0 SV Jss 


dS x A 


Laplacian operator 


V 2 = v.v = 


d 2 d 2 d 2 
dx 2 dy 2 dz 2 


Vector identities 

v(0V’) = 

V.(0A) = A.V</> + 0V. A 

V x (0A) = V0 x A + 0V x A 
V.(A x B) = (V x A).B - A.(V x B) 


V(0 + 0) = V0 + V0 
V.(A+B) = V. A + V.B 
V x (A+B) = VxA + VxB 


V x (V0) = 0 V.(V x A) = 0 

V x (V x A) = V(V.A) - V 2 A 


v 2 0 + if>) = v 2 </> + v 2 0 v 2 OVO = 0V 2 v> + 2 v^.vV’ + V’V 2 ^ 

V(A.B) = (A.V)B + (B.V)A + A x (V x B) + B x (V x A) 

V x (A x B) = A(V.B) - B(V.A) + (B.V)A - (A.V)B 


Cylindrical polar coordinates 


V0 = 

v 2 v> 


dR R d(p dz 
1 d . dip . 1 < 9 2 0 

' RdR^dR’ + Itflty 2 


+ 


d 2 ip 

dz 2 


V.A 


]__d_ 

RdR 


(RA r ) + 


i dAtp 

R d(j) 


+ 


dA z 

dz 


V x A 


( 1 dA z 
[ R dcj) 



dAn 

dz 




9 A r 

d(j> 


k 
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Spherical polar coordinates 


0 


„ , , 1 : 1 5 -ip 

V0 — T~ r H -57-0 H : 7777 

ar r at) r sm d ocp 

VV = — — (r 2 ^^) + 1 d(dn^) + 

r 2 dr dr r 2 sin 9 89 " 86 r * sm 


1 


d 2 ip 


2 ouA ^ Qfjfi 


. 1 d , 9 . , 1 d , . „ . . 

V.A = — — (r 2 A r ) + — - — (sm6»A e ) + 
r z ar r sm t/ at/ 


1 dA^ 

rsinh d0 


V x A = 


r sm t 




r + 


1 


r sin 9 


8A 


8 


90 - sin V' 4 * ) 


+ - 


d dA r 

aF (r/W “ -m 


4> 


General orthogonal curvilinear coordinates 

Transformation from x = (x\,X2,x 3 , ) to u = (111,112,113) 
where xi = Xi(m, U2, u 3 ) i = 1 , 2,3 

dx = (duihiei,du2h2e,du3Ji3e3) 
where hi are scale factors (metric coefficients) 
e,- are unit vectors 


dx 


duj 


i = 1 , 2,3 


and hidi = hi = 

UUi 

Element of arc dl = [h 2 (dm) 2 + h 2 (du 2 ) 2 + h 2 (du 3 ) 2 ~\ 
Area element on ui-coordinate surface dA = h 2 h 3 du 2 du 3 
Volume element dV = h\h2h3du\du2du3 


Vector calculus 


1 dd „ 1 1 dtp 

grad'd = V0 = — 77— ei + — 77— e 2 + —77—63 


hi dm li2 du2 

A = Aiei + A2e2 + A3e3 


h:i du 3 


divA = V.A = 


1 


h:\h2 h.i 




curlA = V x A = 



hie 1 

/1262 

h 3 e 3 

1 

d 

d 

8 

hih 2 h 3 

du 1 

8 u 2 

8 u 3 


/11A1 

h 2 A 2 

h 3 A 3 


v 2 v> = 


1 

' 8 

h 2 h 3 dip 

8 

_| 

h 3 h\ dil’ 

8 

_| 

h 1 h,2 dip 

hi h 2 h 3 

8u 1 

hi dm 

8U2 

h 2 du 2 _ 

du 3 

h 3 du 3 \ 
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7.2 Integral theorems of the vector calculus 


V is the volume bounded by the closed surface S 

Gauss’s divergence theorem 


divF dV = / F.dS 


tv 


Green’s theorem 

First form: 


'V 


+ grad^>.grad;/’)dF = / (/’gradt/’.dS 


Second form: 


j W 2 '!’ ~ i’V 2 <j))dv = J(<f > ^ - *P^)dS 

C is the curve bounding the open surface S 


Stokes’s theorem 


curlF.dS 




F.dr 
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Chapter 8 


Tables of Transforms 

8.1 The Laplace transform 


mm 

coo 

= / e~ st f(t)dt or 

COO 

C[x(t)\ = 


Jo 

Jo 


= F(s) 

= x(s ) 

m 

= C-^Fia)] 

o 

II 

43 s 

e 


Table of Laplace transforms 1 


m c[f{t)} 

C~ l [F{s)\ F(s ) 


1 

1 

s 

(8.1) 

t 

1 

(8.2) 

t n 

n! 

s n+i 

(8.3) 

(n a positive integer) 



t x { A > 1) 

r (A + 1) 

s A+i 

(8.4) 

e -ai 

1 

s + a 

(8.5) 

sin bt 

b 

s 2 + b 2 

(8.6) 

cos bt 

s 

s 2 + b 2 

(8.7) 

sin (at + 4>) 

s sin 4> + a cos 4> 

(8.8) 

s 2 + a 2 



continued. . . 
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Table of Laplace transforms 2 


m 

sinh bt 
cosh bt 
t sin bt 

t cos bt 
sin bt — bt cos bt 
sin bt + bt cos bt 
e~ at t n 

( n a positive integer) 
e~ at sin bt 

e~ at cos bt 

e~ at sinh bt 

e~ at cosh bt 

H(t ) 
H(t — a) 
H(t) - H(t - a) 

m 

5(t — a) 

[t] 

sin at 
t 


tmt)} 

F(s) 

b_ 

s 2 — b 2 
s 

s 2 — b 2 
2 bs 

( s 2 + b 2 ) 2 
s 2 - b 2 
(s 2 + b 2 ) 2 
2 b 3 

(s 2 + b 2 ) 2 
2 bs 2 

(s 2 + 6 2 ) 2 

n! 

(s + a) n+1 


b 

(s + a) 2 + b 2 
s + a 

(. s + a) 2 + b 2 
b 

(. s + a) 2 — b 2 
s + a 

( s + a) 2 — b 2 

1 

s 


s 

1 - e~ as 

s 

1 

e~ as 


s(l - e s ) 

-i a 

tan — 

s 


(8.9) 

(8.10) 

( 8 . 11 ) 

( 8 . 12 ) 

(8.13) 

(8.14) 

(8.15) 

(8.16) 

(8.17) 

(8.18) 

(8.19) 

( 8 . 20 ) 

( 8 . 21 ) 

( 8 . 22 ) 

(8.23) 

(8.24) 

(8.25) 

(8.26) 
continued. . . 
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Table of Laplace transforms 3 


m 

C^lFis)] 


F(s) 


cos 


2 fat 


erfc 


\frri 

a 

2 ft 


“ e -“ 2 /4t 


e ~ a / S 

~V*~ 

s —ay/s 

S 


2\J 1 Tt 3 


(8.27) 

(8.28) 
(8.29) 


g -a 2 /4t 

e ~a^s 

(8.30) 

ffi 

f~S 

a2 ! 4t aerfc ( “ ' ) 

g CL'^S 

(8.31) 

\2 ft) 

ff 

2(1 -a 2 t 2 

f^ 

erfc ( — 'j e s2 t 4a2 

\2 a) 

(8.32) 

a f(t) + (3g{t) 

aF(s) + (3G(s) 

(8.33) 

'e' 

1 

-to 

's' 

1 

-to 

e~ as F(s) 

(8.34) 

e~ at f(t) 

F(s + a) 

(8.35) 

fit ) 

sF(s) - f( 0) 

(8.36) 

fit) 

s 2 F(s) - sf( 0) - f( 0) 

(8.37) 

f (n \t) 

s n F(s) - s n “ 1 /( 0) - s n ~ 2 f{ 0) 

- 


o' 

e 

‘+-» 

1 

(8.38) 

f f {u)du 
Jo 

F(s) 

s 

(8.39) 

-tf{t) 

F'(s) 

(8.40) 

i-trm 

F( n f s ) 

(8.41) 

no 

t. 

roo 

J F(u)du 

(8.42) 

[ f(u)g(t - u)du 

F(s)G(s) 

(8.43) 


(convolution integral) 
f(t + T) = f(t) 

g ofrf Mabt) 

P polynomial of degree less than n 


1 

1 — e~ Ts 


e~ su f(u)du (8.44) 
P(s)/Q(s) (8.45) 


Q(s ) = (s - ai)(s -a 2 )...(s- a n ) 
where ai,a 2 , ■ ■ ■ ,a n are all distinct 
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Limiting theorems 


lim.s^oc F(s) = 0, 


lim^oc f(t) = lim.5^0 sF(s) 


Inversion formula 


1 

2vr j 


n+joo 


e st F(s)ds 


f'y—joc 


f(t) t > 0 
0 t < 0 


Stehfest numerical inversion 

Given f(s), the Laplace transform of f(t), seek the value f(T) for a specific 
value t = T. 

In 2 

Choose Sj = j— j = 1, 2, . . . , M where M is even. 

The approximate numerical inversion is given by 


l 9 M 

f{T) « ^J2wjf( Sj ) 

3 = 1 

where the weights, Wj, are given by 


w 


= (" 1 ) 


%+J 


E 

fc=[|(i+T] 


M 
. 2 


kf{2k)\ 

k)\k\{k - l)!(j - k)\(2k-j)\ 


Stehfest’s weights for M = 6,8,10,12 and 14 


M = 6 

M = 8 

M = 10 

M = 12 

M = 14 

1 

-1/3 

1/12 

-1/60 

1/360 

-49 

145/3 

-385/12 

961/60 

-461/72 

366 

-906 

1279 

-1247 

18481/20 

-858 

16394/3 

-46871/3 

82663/3 

-6227627/180 

810 

-43130/3 

505465/6 

-1579685/6 

4862890/9 

-270 

18730 

-473915/2 

13241387/10 

-131950391/30 


-35840/3 

1127735/3 

-58375583/15 

189788326/9 


8960/3 

-1020215/3 

21159859/3 

-2877521087/45 



328125/2 

-16010673/2 

2551951591/20 



-65625/2 

11105661/2 

-2041646257/12 




-10777536/5 

4509824011/30 




1796256/5 

-169184323/2 

824366543/30 

-117766649/30 
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8.2 The Z-transform 


Z [x{n)\ = ^ x(n)z n = x(z), x(n) = Z 1 [x{z)\ 

n — 0 

Table of Z-transforms 

The variable t is related to the sampling interval, T, by t = nT. 


x{ri) 
Z~ l [x(z)\ 


Z [x(n)\ 
x(z) 


S(n) 
1 or H(n) 



n 

n 2 

n 
k 

a'" 

5(n — m) 
e~ an 

sin(an) 

ne~ an 

cos (an) 

ax\ (n) + f3x2{n ) 
x(n — m) 
e~ an x{n) 
a~ n x(n ) 
x(n + 1) 


1 

z 


z - 1 

Tz 

T 2 z(z + 1) 

(* - l ) 3 

z 

{z - l) fc +! 
z 


z — a 

~—m 


z — e aT 
zsm(aT) 


z 2 — 2z cos (aT) + 1 
Tze~ aT 


( z — e ~ aT ) 2 
z 2 — z cos(aT) 

2 2 — 22 cos (aT) + 1 
0 x 1 ( 2 ) + f3x 2 {z) 
z~ m x(z) 
x(ze aT ) 
x(az) 

z [x{z) — x(0)] 


x(n + 2) 2 2 [^( 2 ) — x(l )2 1 — x(0)] 

x(z)y(z) 


n — m) 


22 x ( m )y( 

m = 0 

Limiting theorems 

z — 1 

lim, woo x(n) = lirn^i x(z), lim x(n) = lim x(z) 

Z n— >0 z—>oo 


(8.46) 

(8.47) 

(8.48) 

(8.49) 

(8.50) 

(8.51) 

(8.52) 

(8.53) 

(8.54) 

(8.55) 

(8.56) 

(8.57) 

(8.58) 

(8.59) 

(8.60) 
(8.61) 
(8.62) 

(8.63) 


45 




8.3 The Fourier transform 


F\f(x)} =-= 

V 2vr J — c 

= m 


f(x)e~ 


^‘ |F(?)] = 7sL 

= fix) 


F{£)e jix dt 


Table of Fourier transforms 


fix) 

^ [/(*)] 

F- 1 1 *’®] 

m 

1 

i 

x 

1 

W\ 

y/2 ne~<H{£) 

a — jx 

H{ \x\) - H(\x\ — a) 

[2 sin a £ 

Vvr £ 

e ~x 2 /a 

#" 2/4 

fix) 


fix)) 

-ep(o 


8.4 Fourier sine and cosine transforms 


F S [fix)} 

=fl) 
= ^,(0 

coo 

1 f(x)s'm^xdx 

0 

Fc [/(*)] 

1 <M 1 ££ 

II II 

COO 

f f(x)cos^xdx 

O 

K 1 [^(01 

■Jh 

= fix) 

COO 

1 F s (£) sill £x(i£ 

0 

•^T 1 [^(01 

= fix) 

COO 

1 F c (£) cos £xd£ 

0 


46 



Table of Fourier sine and cosine transforms 


/ 0) 

m) 

/(*) 

^c(e) 

p-X 

[2 1 

p~X 


xe x2 i 1 2 

V 7T (1 + £ 2 ) 

e -e/2 

-X 2 

e x 

V VT (1 + £ 2 ) 

-f 2 

e ? 

sin x 

X 

1 , !+<£ 

, in 

i — <e 

H(x) — H(x — a) 

|2" sin a£ 

V vr £ 

/'(*) 


f'{x) 

«*■(«) - f-m 

/"(*) 

o) 

A*) 

-em) - \Jlm 


8.5 Some periodic functions: Laplace transforms 
and Fourier series 

Graphs of the Laplace transform, F(s), and the Fourier series, F s (x),F c (x) 
or F(x): 

For the Laplace transform it is assumed that fix) = 0 for x < 0. 


1. Square wave (odd function) height h, period 2 a 


j —h —a < x < 0 
( h 0 < x < a 


F(s) = 



4/ly, 1 

7 T ^ ( 2 n + 1 ) 

n=0 v ' 


sin 


(2 n + \)ttx 
a 
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2. Square wave (even function), period 2 a 



—a < x < —a/2 
-a/2 < x < a/2 
a/2 < x < a 


F(s) = — e aS//2 tanh (Jip) 


F c (x) 


(~l) n 

7 r /2n + 1) 

71=0 V ^ 


COS 


(2 n + l)7rx 
a 



3. Triangular wave, period 2 a 


( —hx/a —a < x < 0 
( hx/a 0 < x < a 



4. Sawtooth wave, period 2 a 

f(x) = hx/2a 0 < x < 2a 
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5. Rectified sine wave, period a 


f(x)= hsm(x/a) 0 <x<a 



6. Half-rectified sine wave, period a 


/(®) 


hsm(2x/a) 0 < x < a/2 
0 a/2 < x < a 
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Chapter 9 

Numerical Methods 


9.1 Interpolation 


Polynomial interpolation: Suppose that / is defined on the interval — 1 < x < 
with n + 1 continuous derivatives. If the n-degree polynomial approximation 
at the n +1 points xq , X \ , . . . , x n is given by p n ( x ) then the error at any point 
x £ (—1, 1) is given by 


ft \ t \ ( X ~ X 0 )(x-Xi)... (x-X n ) ,(„+!), f . ,v 

f{x)-p n {x) = 7 — — tt / l + +6 for some ^ G (-1, 1) 

[n + 1 )! 


Lagrange interpolation polynomials 


(■ X - x 0 )(x - Xl) . . . (x - Xj-l)(x - Xj+l) ,..(x-x n ) 
(Xi - Xo )(Xi - Xl) ... (Xi - Xi-l)(Xi - Xj+l) . . . (Xi - x n ) 


The n th degree collocation polynomial through the points 
(xo,yo), (xi,yi) ■ ■ • (x n ,y n ) is given by 


n 

P n {x) = ^ Li(x)yi 

i= 0 


Cubic splines 

The cubic spline interpolating function, S(x), through the points 

(xo, Vo), (x'i , yi) ■ ■ ■ (x n , y n ) has the properties S(xi) = y t i.e. S is continuous 

at (xj, yi). 

S'(x) and S"{x) are continuous. 

S(x) is a cubic polynomial, Sj(x), in each interval [x+Xj+i]. 

y i+ i{x-Xi) yi(x-x i+ i) 

Si(x) = - 


hi 

Mi 

6 


(x - x i+ iY 
hi 

where M,- are found from 


hi 

- hi (x - Xj+l) 


Mi 


i + 1 


x — Xj 


- hi(x - Xi) 
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a 0 o doi 0 0 • • • 0 M 0 

/i 0 2(/io + /ii) h i 0 • • • 0 Mi 

0 hi 2(/ii + h‘2 ) h‘2 ... 0 M 2 

0 • • • 0 h n — 2 2(/i n _2 + h n ~ 1 ) h n — 1 M n _i 

0 • • • 0 0 CL nn — 1 Ann _ _ -M n 

bo 

V 2 -yi _ yi—yo 
hi ho 

V3-y2 _ V2—yi 
h.2 hi 


where hi = Xj+i — Xj. 

Natural spline, Mo = M n = 0 so that 

®00 = G-nn = 1) ®01 = ®nn- 1 = bo = b n = 0 

If the derivative of the underlying function is known at xo and yo then 
a 0 o = 2/r 0 , a 0 i = h 0 , 6 0 = (yi - y 0 )//i 0 - /'(x 0 ) 

®rm- 1 — h n — i, CL nn — 2h n —\, b n — f (X n ) (Un Vn— l)/hn— 1 

9.2 Finite difference operators 

Forward Ay* = y i+1 - y.j Backward Vy* = y* - yi-i 

Central <5y i+ i = y l+ \ - yi Shift Ey % = y i+l 

Average yy i+ i = (y, + y i+ i)/2 Differential Dy t = (dy/dx) Xi 


Relationship between the operators 



E 

A 

8,ii 

V 

hD 

E 

E 

1 + A 

1 + y>8 + e> 2 /2 

(1-V)- 1 

e hD 

A 

E- 1 

A 

+ 5 2 /2 

V(1 - V)- 1 

e hD - 1 

8 

E 2 — E~ 2 

A(1 + A)5 

<5 

V(1 - V)"5 

2 sinh(FF?/2) 

V 

1 - FT 1 

A(l + A)" 1 

fiS - 8 2 / 2 

V 

1 - e~ hD 

hD 

in A 

ln(l + A) 

2 sinh _1 (<5/2) 

- ln(l - V) 

hD 

h 

(F?^ + E~^)/2 


(1 + <5 2 /4) i 


cosh(AD/2) 
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9.3 Non-linear algebraic equations 

Formula iteration method for x = F(x). 

An iterative scheme has the form x n +i = F(x n ), which converges to the root 
x = a, only if | F'(a)\ < 1. The convergence is hrst order; i.e. if |F r ' / (ic)| < M 
for x near a, then the error at each stage satisfies |e n +i| < M|e n |. 

Newton- Raphson method for f(x) = 0. 


•En+1 — %n 


f{Xn) 

f'(x n ) 


The method converges to the root x = a if xq is sufficiently close to a. The 
convergence is second order i.e. if \ f"(x)\ < 2 M for x near a then the error 
at each stage satisfies |e n+ i| < M\e n \ 2 . 


Newton’s method for a system of equations f(x) = 0, where 
f = {/i(x),/ 2 (x), . . • ,/ n (x)} and x = {xi,x 2 ,.. -,x n }. 

x (n+i) = x (n) _ 

where J n is the Jacobian matrix given by 


(J 


n )ij 


' dfj 
Pxj 



9.4 Numerical integration 

The region of integration is subdivided into n equal intervals of width 
h = (b — a) /n, Xi = a + ih and /* = f(xi). 

Trapezoidal rule 

f{x)dx « h{(fo + fn )/2 + (/l + /2 + • • • + /re— l)} 

The error, e n , is bounded as follows: 

h 2 

|e n | < (b — a) — Mi, where Mi = max |/ // (.x)| 

12 x£[a,b] 

Simpson’s rule (n must be even) 

J f{x)dx ~ — {/o + fn + 4(/i + fz + • • • + fn- 1) + 2(/2 + fi + • • • + fn- 2 )} 
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The error, e n , is bounded as follows: 

|e n | <(b- a)— -M 2 , where M 2 = max |/ (lv) (x)| 

loU x€[a,b\ 

Gauss quadrature 

Whenever f(x) is known explicitly Gauss quadrature can be employed: 


pi n 

/ f(x)dx « 

J ~ l i = 1 


The error, e n , 

M 


is bounded as follows 
2 2 n +1 ( n !) 4 


< 


(2n + l)[(2?r)!] 3 


M 3 where M 3 


max |/( 2 n) (.x)| 
xe[-i,i] 


Gauss quadrature: integration points and weights 


n 


±Xi 



UJi 


1 

0.00000 

00000 

00000 

2.00000 

00000 

00000 

2 

0.57735 

02691 

89626 

1.00000 

00000 

00000 

3 

0.00000 

00000 

00000 

0.88888 

88888 

88889 


0.77459 

66692 

41483 

0.55555 

55555 

55556 

4 

0.33998 

10435 

84856 

0.65214 

51548 

62546 


0.86113 

63115 

94053 

0.34785 

48451 

37454 

5 

0.00000 

00000 

00000 

0.56888 

88888 

88889 


0.53846 

93101 

05683 

0.47862 

86704 

99366 


0.90617 

98459 

38664 

0.23692 

68850 

56189 

6 

0.23861 

91860 

83197 

0.46791 

39345 

72691 


0.66120 

93864 

66265 

0.36076 

15730 

48139 


0.93246 

95142 

03152 

0.17132 

44923 

79170 


The Xi are the positive zeros of the Legendre polynomial P n (x). The 
formula integrates exactly all polynomials of degree 2 n — 1. Any interval 
a < X < b can be transformed to the interval — 1 < x < 1 by the change of 
variable 

X = (a + b)/2 + (b — a)x / 2 

For double integrals: 


/ (x, y)dxdy 


1-1 J - 1 


n n 

EE 

i = 1 j = 1 
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9.5 Numerical solution of ordinary differential equations 

First order equations 

y' = f(x,y) with y(x 0 ) = 2/0 

The equation is solved successively at discrete points x\, X2, x 3 , . . . 

where Xj+i — Xi = h. 

y n is the exact value y(x n ). 

Y n is the approximation to y n , used in the recurrence relation. 

Single-step methods 

Euler’s method 

Yi+^Yi + hfix^Yi)-, Y 0 = y 0 
h 2 

The error is —y"{fi) where x* < £ < Xi + \. 

Modified Euler method 

Y i+ i =Yi + \ Yi ) + fixi+uY^)} , Yq = 2/ 0 

where Y^_ 1 = Yi + hf{xi , ¥>) 

h 3 

The error is ——y'"[fi) where xi < £ < Xi + \. 

Runge-Kutta (fourth order formulae ) 

Yi + 1 = Yi + — [k\ + 2 k 2 + 2 k 3 + kfi), Yq = y 3 
o 


where ki = hf(x il Yi), k 2 = hf(xi + h/2, Yi + Aq/2), 

h = hf (xi + h/2, Yi + k 2 / 2), fc 4 = hf(xi + h,Yi + k 3 ). 

The actual form of the local truncation error, ej+ 1 , is extremely complicated; 
it is, however, of the order hfi . 


Multi-step methods (Predictor-corrector) 
Milne-Simpson method 

A h 

Predict Y^ = Y t _ :i + —[2/* - /*_ 1 + 2 /,_ 2 ] 

Correct = Y t _ 1 + ^ [/*— 1 + 4/j + /^r] 

The error in the predictor is — /i 5 2/ v ^(£i) 

h 5 

and the error in the corrector is — gg y^(&) 
where Xi- 3 < £1 < x i+1 , Xi - 1 < 6 < ®i+i- 
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Adams-Moulton method 

Predict Y ) ^ ] = Yj, + ^-[55 /, - 59/*_i + 37/j_ 2 - 9/j_ 3 ] 

Correct = y + 7 ^[fi -2 ~ 5f[-i + 19 fi + 9/j+i] 

The error in the predictor is —h 5 y^ v \^ i) 

h 5 

and the error in the corrector is — — y^ v ^(^ 2 ) 
where Xi < £1 < Xi + \ and Xi < £2 < ®i+ 1 - 


Higher order equations 

The methods given above are easily generalised for use with a system of 
first order equations, thus we consider only the first order equation in one 
dependent variable. 

The n th order equation y ^ = f(x,y,y', . . . may be reduced to 

the system of first order equations 

2/1=2/, 2/2 = 2/1, 2/3 = 2 /!, ••• Vn—l = y'n—2 

y'n = /(z, 2/1, 2/2,. ■ ■ , 2 /n) 

i.e. y' = f(x,y) 

The corresponding Euler recurrence relation is 

Yj+i = Yj + yj, Y 0 = y 0 


9.6 Systems of linear equations, n x n 

011^1+ O 12 X 2 T • • • T 0\ nXn — h\ 

U21*Xl+ 0,22X2 + . . . + 02 n X n = h'2 

in matrix form Ax=b 

^niXiT o n 2 X 2 T • • • T o nn x n — h n 

Direct method 

Gauss elimination , with partial pivoting 
Computational procedure, (k = 1 . . . (n — 1)): 

1. Rearrange the equations so that 
\akk\ > Wik\, i = (k + l)(l)n 

2. Compute 

hik O'ik / ^kki i {k T 1) ••• n 

3. Compute 

a'ij = aij + ytkOkj, j = (k + 1) . . . n 1 
hi = hi + Hikbk J 
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4. Compute 

x n = b' n /a‘, 


i /„/ 
nn 


Xj = 


h 'i ~ a ij X i ) / a 'ii ’ * = ( n - X ) ' ' ' 1 

J=*+l 


Indirect methods 

Jacobi 

An iterative scheme is given by 


x. 


(r+l) _ 1 


j = i 


' j 


(r) 


* = 1, 2, . . . , n 


In matrix form 


X ( r+1 ) = b - [L + U]x< r ) 


Gauss-Seidel 

An iterative scheme is given by 


x. 


(r+l) _ 1 


bi-^2 a^ +1) 

3 = 1 


22 a ij x 2 ) ’ * = 

j=i + 1 


A sufficient condition for convergence of both methods is that the matrix A 

is diagonally dominant 

i.e. 




'y ^ I a ij 1 1 * — 1,2 , ,n 


j = i 

3 + 1 


In matrix form 


where 


x 


(r+i) = b _ Lx (r+a) _ UX (r) 


0 

«21 
a 22 
031 
033 


0 

0 

032 

033 


gnl Qn2 
(Jnn &nn 


0 

0 

0 


0 

0 

0 

®nn — 1 


0 

0 

0 


,U = 


0 21a 

Oil 

0 0 


ai3 


(Jin 

an 

CL23 


an 

a 2n 

CL22 


CL22 

0 

0 

&n — In 

(Jn — In — 

0 

0 

0 


and b = [bi/an 62/022 • • • b n /a r 
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The Gauss-Seidel process converges if and only if all the eigenvalues of 
the matrix [I + L]-!U have modulus less than one. 

Successive over-relaxation (SOR) 

The SOR iterative scheme is 

x (r+l) = x (r) + w ( b _ Lx (r+1) _ X W _ Ux M) i 

where 1 < to < 2 for over-relaxation, and u = 1 for Gauss-Seidel. 

The SOR process converges if and only if all the eigenvalues of the matrix 
[I + o;L] _1 [(l — cu)I — ojU] have modulus less than one. 

9.7 Chebyshev polynomials 

Recurrence relation 

T n+1 (x) = 2 xT n (x) - T n _i(x); T 0 (x) = 1, Ti(x) = x 
T n (x) = cos(ncos _1 x) (coefficient of x n is 2 n ~ l ) 

T^x) = 2x 2 — 1, Ts(x) = 4x 3 — 3x, T^x) = 8.x 4 — 8x 2 + 1, 

T 5 (x) = 16x 5 — 20x 3 + 5x 

The polynomial (l/2 n ^ 1 )T„(x) has a smaller upper bound to its magnitude 
over [—1,1] than that of any other polynomial with leading term x n . 


9.8 Numerical eigenvalues and eigenvectors 

The matrix eigenvalue problem is given by, see page 13, 

Ax = Ax 


Power method 

(for a non-repeated dominant eigenvalue Ai, and eigenvector xi). 

Let zq be an arbitrary vector, but not xi,X2, . . . ,x n , then the iterative pro- 
cess Zj+i = r-Az i (where fc* is the element with the largest absolute value 
in z j) converges, with ki — > Ai and z t — > xi, provided that A has n linearly 
independent eigenvectors. The convergence will be slow if | A 2 1 ~ |Ai|. 
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Inverse iteration 

(for the eigenvalue A closest to p and eigenvector x) 

Let zq be an arbitrary vector, then the iterative process 

Zi+i = -7 [A -pI] _1 Zj 

Ki 

converges with fc* — * and z * — > x, where ki is the element with the 

largest absolute value in z,. 

Jacobi method 

(for symmetric matrices) 

Suppose that A is diagonalised by using a sequence of orthogonal transfor- 
mations 


D = . . . T£T£ ATiT 2 . . . T k _iT k = M* AM, say, 

then the columns of M are the eigenvectors and the diagonal of D comprises 
the corresponding eigenvalues. 

Computational procedure: 

1. Locate largest off-diagonal element a pq , say. 

2. Compute 9, where tan 26 = 2 a pq /(a qq — a pp ), \9\ < 7r/4. 

3. Compute new elements in rows p and q 

Q'pp — (tan 9 )u pq (i qq — ciqq H - (tan 9)a pq 

CLpq 0 

d'pj = (cos 9)a p j — (sin0)a g:; ' a! qj = (sin0)a P:; ' + (cos 6)a q j 

4. If any off-diagonal element is non-zero, return to 1. 

LR method 

(for all the eigenvalues) 

Form the sequence Ao = A, A r+ i = U r L r , where A r = L r U r with 
L r lower triangular, all diagonal elements are equal to 1 and U r is upper 
triangular. 

For suitable A the sequence converges to an upper triangular matrix whose 
diagonal elements are the eigenvalues of A, arranged in order of decreasing 
modulus. 

If the method fails to converge try applying it again to A + pi where p is a 
suitable real number. Convergence in this case will be to A* + p. 
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9.9 Least squares approximation 

Given the data set {(xo, Do), (xi, yi ), . . . , (xjv, Un)}, the least squares best 
fit n th degree polynomial (n < N — 1) is 

y = a 0 + ci\x + a 2 x 2 + . . . + a n x n 

where the coefficients a* are found from the normal equations 


’ n Y.*i Ex? ••• Ex” 


Oo 


E Vi 

Ex. Ex? Ex? ■■■ E x? +1 


ai 


E x iVi 

E 4 ■■■ E x? +2 


a 2 

= 

YAvi 

symmetric : 





E x l n . 




. E x iVi . 
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Chapter 10 

Statistics 


10.1 Sample statistics 


The following definitions are concerned with the data set xi,x 2 , ■ • • ,x n . 
Where appropriate, the data value x t occurs with frequency 


Sample Mean 

(arithmetic) 
or Average 


Population variance 


Standard deviation 


Sample variance 


Pooled variance 
Standard error of mean 


x = 


a 2 = 


-t n k k 

~H xi or 2-WZ^ 


i = 1 
n 


i = 1 


i= 1 


Tl 


B 

i= 1 


x,- — x) 2 = 


E*? - 


2 (E 


n 


n 


or 


a 2 = 


E* 

2=1 


■E« 


f 2 __ (E fi x i) 2 

J iXi r 

-\2 _ B Ji 


Xj — X ) = 


7 = 1 


(7 


•S 2 = 


1 " 

^tE( 


n — 1 


x,- — x) 2 = 


Eb - 


E/i 


2 (E^ 2 


n 


i=l 


(n- !) 


2 _ (ni - l)sf + (n 2 - l)s 2 
Sp ni + n 2 - 2 

<7 


n 
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10.2 Regression and correlation 


Pearson’s Product-moment 
Correlation Coefficient of 
n sample pairs (x t ,yi) 


r = 


E( x i — x) (jjj -y) 

VIE (Xi-xyum-y) 2 ] 
E 

n 


— 


(Ew) 


21 


n 


n 


Least squares estimates 
a and b in the fitted 
regression line y = a + bx 


. = E( x i - x )(yi - y ) 

E(Xi-x) 2 
a = y — bx 


E XiVi - 


Ex 


(Exj)(Eyi) 

n 

iExi ) 2 


Residual Variance 
in simple linear regression 


^2 ~ a ~ bx rf 


SSE 
n — 2 


Spearman’s rank correlation 
coefficient of n pairs 
of sample rankings ( Xi,yi ) 


r s = 1 - 


n(n 2 — 1) 


where di = X{ — yi 


10.3 Distributions 


Discrete uniform distribution 

k: the possible number of values 


Probability function 


Mean 


- i 

X — X\, X 2 5 • • • 5 %k 

k + 1 
2 


Variance 


k 2 — 1 
12 


Hypergeometric distribution 

V: population size 
n: sample size 

k: number of items in population labelled “success” 
Probability function P[x) = 
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Mean 


Variance 


nk 

V 

k 



( N-n\ 

\N-l) 


Binomial distribution 

n: number of trials 

p: probability of “success” at each trial 
Probability function P[x) = 



p x (l -p) n ~ x 


x = 0, 1, . . . , n 


Mean np 

Variance np(l — p) 


Negative binomial distribution 

r: no. of successes required 

x : no. of the trial when the r th “success” occurs 

or alternatively, if c is the number of “failures” preceding the r th 

Probability function P{x) = f ^ ^ \p r {l—p) x ~ r 


Mean 

Variance 

Probability function 
the general term of the 
Mean 

Variance 


x = r,r + 1, . . . , oo 
r 

P 

r(l — p) 

p2 

P(c)= ^ C + r c ~ 1 ^ p r q c 

binomial expansions of p r (l — q)~ r 
r(l — p) 

P 

r(l — p) 


Geometric distribution 

p: probabililty of “success” at each trial 
Probability function P(x) = p( 1 — p) x ~ l 


Mean 


re = 1 , 2 ,... 
1 

P 


Variance 


1 — p 


p- 


“success” 
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Poisson distribution 

m: average number of random events in a given interval 

g TYlyyyX 

Probability function P(x) = — 

' 7*1 

x = 0, 1, . . 

Mean m 

Variance m 


Continuous uniform distribution 

a: minimum value 
(3: maximum value 

Probability function f(x;a,/3) = — — - 

a < x < /3 


Mean 


Variance 


f3 -\~ OL 
2 

(P - q ) 2 

12 


Exponential distribution 

m: average number of random events in a given interval 
Probability function f{x\m ) = me~ mx x > 0 

Mean — 

m 


Variance 


1 

m 2 


Normal distribution 

/j: mean 
cr 2 : variance 

Probability function 


Mean 

Variance 


f(x-fi,a 2 ) = — /== exp 
(iv27r 

— oo < x < oo 



V b 
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Chapter 11 

S. I. Units (Systeme 
International d’Unites) 

11.1 Fundamental units 

The kilogram (kg) is defined as the mass of the international prototype 
platinum- iridium cylinder kept at Sevres. 

The second (s) is defined as the time taken by 9 192 631 770 periods of the 
radiation from the transition between the two hyperfine levels of the ground 
state of the atom Caesium- 133. 

The metre (m) is defined as the path length of light in a vacuum during an 
instant of 1/(2.99792458 xlO 8 ) seconds. 

The mole (mol) is defined as the amount of substance of a system which 
contains as many elementary entities as there are atoms in 0.012 kilograms 
of the isotope Carbon-12. 

The ampere (A) is defined as that constant current which, if maintained in 
each of two infinitely long straight parallel wires of negligible cross-section 
placed 1 metre apart in vacuo, would produce, between the wires, a force 
2 X 10 -7 newtons per metre length. 

The kelvin (K) is defined as the fraction 1/273.16 of the thermodynamic 
temperature of the triple point of water. 

The candela (cd) is defined as the luminous intensity from a source of 
monochromatic radiation of frequency 540xl0 12 Hz which has a radiant in- 
tensity of 1/683 Watts per steradian in a given direction. 
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11.2 S. I. Prefixes and multiplication factors 


Factor 

Prefix 

Symbol 

10 24 

yotta 

Y 

10 21 

zetta 

Z 

10 18 

exa 

E 

10 15 

peta 

P 

10 12 

tera 

T 

10 9 

giga 

G 

10 6 

mega 

M 

10 3 

kilo 

k 

10 2 

hecto 

h 

10 1 

deca 

da 

KT 1 

deci 

d 

10 -2 

centi 

c 

10“ 3 

milli 

m 

10 -6 

micro 

M 

1(T 9 

nano 

n 

1(T 12 

pico 

P 

10“ 15 

femto 

f 

10- 18 

atto 

a 

10- 21 

zepto 

z 

10- 24 

yocto 

y 


11.3 Basic and derived units 


Table of basic units 


Physical quantity 

Dimensions 

S. I. unit Symbol 

mass 

M 

kilogram kg 

length 

L 

metre 

m 

time 

T 

second 

s 

amount of substance dimensionless mole 

mol 

electric current 

I 

ampere 

A 

temperature 

6 

kelvin 

K 

luminous intensity 

C 

candela 

cd 


Table of supplimentary units 

Physical quantity 

Dimensions 

S. I. unit 

Symbol 

angle 

dimensionless 

radian 

rad 



degree 

O 

solid angle 

dimensionless 

steradian 

sr 


2ir rad = 360° ; 

1 rad = 57.296° 
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Table of derived units 


Physical quantity 
area 
volume 
frequency 
speed 

acceleration 
angular speed 
angular acceleration 
density 
momentum 
moment of inertia 
angular momentum 
force 

torque, moment of force 

energy, work 

pressure, stress 

power, radiant flux 

viscosity (dynamic) 

surface tension 

electric charge 

electric potential 

electric resistance 

electric conductance 

electric field strengh 

electric charge density 

electric flux density 

capacitance 

magnetic flux 

magnetic flux density 

inductance 

permittivity(e) 

permeability(/j) 

quantity of heat 

heat capacity, entropy 

specific heat capacity, specific entropy 

radiant flux density 

thermal conductivity 

latent heat 

molar energy 

molar entropy, molar heat capacity 
luminous flux(lcd=ll m sr -1 ) 
illuminance 


S. I. unit (symbol) 


I? 

L 3 
T~ i 

LT~ 1 
LT ~ 2 

T -i 

T~ 2 
ML ~ 3 
MLT - 1 
ML 2 
ML 2 T ~ 1 
MLT ~ 2 

ml 2 t ~ 2 
ml 2 t ~ 2 

ML _i r _ 2 

ml 2 t ~ 3 

ML _1 T _1 
Mr 2 
TI 

ML 2 T~ 3 I ~ 1 
ML 2 T~ 3 I ~ 2 
M-'L~ 2 T 3 I 2 

mlt~ 3 i ~ 1 
l~ 3 ti 
l~ 2 ti 
m~ 1 l~ 2 t 4 i 2 

M L 2 T~ 2 1~ 2 

mt -2 7 -l 

ml 2 t~ 2 i ~ 2 

M-'L~ 3 T 4 I 2 

mlt~ 2 i ~ 2 
ml 2 t ~ 2 

ML 2 T- 2 0~ l 

L 2 T - 20~1 

MT ~ 3 
MLT- 3 6 ~ 1 
L 2 T~ 2 

ml 2 t ~ 2 

ML 2 T~ 2 0- 1 

ml 2 t ~ 3 

MT 1 -' 2 


Dimensions 


m 

m 

hertz (Hz) 

ms -1 

_2 

ms 

rads -1 

rads - 2 

kgm -3 

kgrns -1 

kgm 2 

kgm 2 s _1 

newton(N) 

Nm 

joule(J)Nm 

pascal(Pa)Nm- 2 

watt(W)Js _1 

Pas 

Nm -1 

coulomb (C) 

volt(V)WA _1 

ohm(H)VA _1 

siemens (S) 

Vm _1 

Cm -3 

Cm- 2 

farad(F)AsV _1 

weber(Wb) 

tesla(T)Wbm- 2 

henry (H)WbA -1 

Fm _1 

Hm _1 

joule(J) 

JR" 1 

Jkg _1 K~ 1 

Wnf 2 

Wm _1 K _1 

Jkg' 1 

Jrnol- 1 

Jmol - 1 !^- 1 

lumen(l m ) 

lux(l x ) 
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11.4 Values of some physical constants 


Quantity 

Symbol 

Value 

speed of light in a vacuum 

o 

II 

/ nT' 

o 

o 

Mi- 

2.99792458 x K^ms” 1 

permeability of a vacuum 

M o 

Air X lCT'Hm -1 

permittivity of a vacuum 

eo 

8.854187817 x 10“ 12 FM' 

elementary charge 

e 

1.60217733 x 10” 19 C 

Coulomb force constant 

k e 

8.9875 x 10 9 Nm 2 C“ 2 

fine structure constant 

a = e 2 / (47reo He) 

1/137.0 

gravitational constant 

G 

6.67258 x 10“ 11 m 3 kg~ 1 s 

atomic mass unit 

u 

1.6605402 x 10 _27 kg 

energy equivalent of u 


931.494MeV 

rest mass of electron 

m e 

9.1093897 x 10“ 31 kg 
0.0005486u 

electron rest energy 

m e c 2 

0.5109991MeV 

rest mass of proton 

m p 

1.6726231 x 10~ 27 kg 
1.007276u 

proton rest energy 

m p c 2 

938.27231MeV 

rest mass of neutron 

m n 

1.6749286 x 10~ 27 kg 
1.008665u 

neutron rest energy 

m n c 2 

939.56563MeV 

Planck’s constant 

h 

6.626075 x 10“ 34 Js 


h = h/ (27 r) 

1.05457266 x 10“ 34 Js 

Rydberg energy 

Roo = a 2 m e c 2 /2 

13.61eV 

Bohr radius 

«o 

5.29177249 x 10“ n m 

Compton wavelength of electron 

Ac 

2.42631058 x 10“ 12 m 

Avogadro constant 

Na 

6.0221367 x lO^mol" 1 

Boltzmann constant 

k 

1.380658 x 10~ 23 JK _1 

Stefan-Boltzmann constant 

a 

5.67052 x 10 _8 Wm _2 K~ 

Wien’s (displacement law) constant 


2.8978 x 10 _3 mK 

Bohr constant 

Mb 

9.2740154 x lO'^JT” 1 

nuclear magneton 

Mn 

5.05057866 x lO-^JT" 1 

triple point temperature 

T t 

273. 16K 

molar gas constant 

R 

S.SISJK -1 !^!^ 1 

micron 

Mm 

10” 6 m 

angstrom 

A 

10” 10 m 

femtometre or fermi 

fm 

10" 15 m 

barn 

11.5 Useful masses 

b 

10 _28 m 2 


The atomic mass unit, u, is defined to be yg of the mass of the Carbon-12 
atom. 
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Table of particle masses 


particle/atom mass 


p to p = 1.6726231 x 10 _27 kg 

= 1.007276u 
= 938.272MeV/c 2 

e m e = 9.1093897 x 10 -31 kg 

= 0.0005486u 
= 0.511MeV/c 2 

n m n = 1.6749286 x 10~ 2 'kg 

= 1.008665u 
= 939.566MeV/c 2 


}H 

1.007825u 

?H 

2.014102u 

?H 

3.016049u 

2 He 

3.016029u 

2 He 

4.002603u 

iLi 

6.015121u 

7 Li 

7.016003u 

5 2 B 

12.014352u 

12p 

6 ^ 

12.000000U 

13p 

6 ^ 

13.003355U 

14 p 

6 ^ 

14.003242U 

42 n 

12.018613ii 

fN 

13.005738ii 

^ 4 N 

14.003074u 

iKr 

91.8973U 

56 1 Ha 

140.9139u 

i°Po 

209.982848U 

IfPo 

218.008965u 

220 p 

86 ™ 

220.011369u 

li 2 Hn 

222.017574u 

| 4 Ra 

224.020187u 

2 | 6 Ra 

226.025402u 

90 ? Th 

227.027701u 

90 8 Th 

228.028716u 

io°Th 

230.033127u 

231 Pa 

91 rd 

231.035880u 

90 4 Th 

234.043593u 

235 it 

92 u 

235.043924u 

238 tt 

92 u 

238.050784u 

l 5 Np 

235.044057u 
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11.6 Astronomical constants 


Mass of Earth 

m E 

5.976 x 10 24 kg 

Radius of Earth (equatorial) 

Re 

6.378 x 10 6 m 

Gravity at Earth’s surface 

9 

9. 80665ms -2 

Mass of Sun 

M 0 

1.989 x 10 30 kg 

Radius of Sun 

Rq 

6.9599 x 10 8 m 

Solar effective temperature 

T e 

5800K 

Luminosity of Sun 

L© 

3.826 x 10 26 W 

Astronomical unit 

AU 

1.496 x 10 11 m 

Parsec 

pc 

3.086 x 10 16 m 

Jansky 

Jy 

10 -26 Wm -2 Hz 

Tropical year 


3.1557 x 10 7 s 

Standard atmosphere 

atm 

101325Pa 


11.7 Mathematical constants 


pi (Archimedes’ constant) 
exponential constant 
Apery’s constant 
Catalan’s constant 
Dottie’s number 
Euler’s constant 
Feigenbaum’s constant 
Feigenbaum’s constant 
Gelfond’s constant 
Gibb’s constant 
Golden mean 
Khintchine’s constant 
omega constant 
parabolic constant 
plastic constant 
Sierpiriski’s constant 
Trott’s constant 
Wallis’s constant 


7 r 

3.141 

592 

653 

793 

e 

2.718 

281 

828 

459 

C(3) 

1.202 

056 

903 

160 

K 

0.915 

965 

594 

177 

r 

0.739 

085 

133 

215 

7 

0.577 

215 

664 

902 

a 

2.502 

907 

875 

096 

6 

4.669 

201 

609 

102 

e n 

23.140 692 632 78 

G 

1.851 

937 

051 

982 

</> 

1.618 

033 

988 

750 

K 

2.685 

452 

001 

065 

n 

0.567 

143 

290 

410 

P2 

2.295 

587 

149 

393 

p 

1.324 

717 

957 

245 

I< 

2.584 

981 

759 

579 

T 2 

0.108 

410 

151 

223 

w 

2.094 

551 

481 

542 
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11.8 The Greek alphabet 


Letter lowercase uppercase 


Alpha 

a 

A 

Beta 

P 

B 

Gamma 

7 

r 

Delta 

S 

A 

Epsilon 

e 

E 

Zeta 

c 

Z 

Eta 

V 

H 

Theta 

9 

0 

Iota 

L 

I 

Kappa 

K 

K 

Lambda 

A 

A 

Mu 

9 

M 

Nu 

V 

N 

Xi 

t 

77 

Omicron 

o 

0 

Pi 

7 r 

n 

Rho 

P 

p 

Sigma 

a 

E 

Tau 

T 

T 

Upsilon 

V 

Y 

Phi 


$ 

Chi 

X 

X 

Psi 

V’ 

T 

Omega 

u>, w 

n 
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Index 


Acceleration vector, 31 
Adams-Moulton method, 55 
Ampere, 64 
Analytic function, 9 
Apery’s constant, 69 
Arc length of a curve, 27 
Archimedes’ constant, 69 
Argand diagram, 8 
Argument, 8 
Arithmetic mean, 9 
Astronomical constants, 69 
Auxiliary equation, 33 
Average, 60 

Bernoulli’s inequality, 9 
Bessel function, 23 
Bessel’s equation, 23 
Beta function, 22 
Binomial series, 25 
Binomial theorem, 3 

Candela, 64 
Catalan’s constant, 69 
Cauchy Schwarz inequality, 9 
Cauchy-Riemann equations, 9 
Centroid, 28 
Chain rule, 17 

ordinary differentiation, 17 
partial differentiation, 17 
Chebyshev polynomials, 22, 57 
Circle, 6, 7 
Circumcircle, 7 
Cofactor, 10 
Combinations, 3 

Complementary error function, 22 
Complementary function, 33 


Complex conjugate, 8 
Complex numbers, 7 
Complex roots, 8 
Complex variable, 9 
Cone, 7 

Conic sections, 6 
Convolution, 43 
Correlation, 61 
Cosine rule, 7 
Cramer’s rule, 11 
Cubic equation, 2 
Cubic splines, 50 
Curl, 37 
Curvature, 23 
Cylinder, 7 

De L’Hopital’s rule, 24 
De Moivre’s theorem, 4, 8 
Degrees and Radians, 3 
Determinant, 10 
Diagonally dominant, 56 
Differentiation, 17 
chain rule, 17 

differentiation of integrals, 18 
Leibniz rule, 17 
linearity rule, 17 
product rule, 17 
quotient rule, 17 
standard forms, 19-21 
total differential, 18 
Directional derivative, 37 
Directrix, 6 
Distributions, 61 
binomial, 62 
continuous uniform, 63 
discrete uniform, 61 
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exponential, 63 
geometric, 62 
hypergeometric, 61 
negative binomial, 62 
normal, 63 
Poisson, 63 
Div, 37 

Dottie’s number, 69 
Dynamics, 31 

Eccentricity, 6 
Eigenvalues, 13 

Eigenvalues and eigenvectors, 13 
inverse iteration, 58 
Jacobi method, 58 
LR method, 58 
power method, 57 
Eigenvectors, 13 
Ellipse, 6, 7 
Error function, 22 
Euler’s constant, 69 
Euler’s formulae, 8 
Euler’s method, 54 
modified, 54 

Exponential constant, 2, 69 
Factorial, 1 

Factors and roots of equations, 2 
Feigenbaum’s constant, 69 
Finite difference operators, 51 
Finite series, 1 
Focus, 6 

Formula iteration, 52 
Fourier series, 35, 47 
complex series, 36 
half-range series, 35 
sine and cosine series, 35 
Fourier transform, 46 

sine and cosine transform, 46 
Fundamental theorem of algebra, 8 
Fundamental units, 64 

Gamma function, 22 
Gauss elimination, 55 
Gauss quadrature, 53 


double integrals, 53 
Gauss’s divergence theorem, 40 
Gauss-Seidel iteration, 56, 57 
Gelfond’s constant, 69 
General conic, 7 

General orthogonal curvilinear coor- 
dinates, 39 
Geometric mean, 9 
Geometry, 5 
Gibb’s constant, 69 
Golden mean, 69 
Grad, 37 

Gram Schmidt process, 12 

Green’s theorem, 40 

Green’s theorem in the plane, 27 

Half-rectified sine wave, 49 
Heron’s formula, 7 
Hessian, 23 
Hyperbola, 6, 7 
Hyperbolic functions, 5 

Imaginary part, 8 
Indices, 1 
Inequalities, 9 
Initial condition, 32 
Initial-value problem, 32, 34 
Inscribed circle, 7 
Integrating factor method, 32 
Integration, 18 
by parts, 18 
by substitution, 18 
definite integrals, 21 
differentiation of integrals, 18 
standard forms, 19-21 
Interpolation, 50 
Intrinsic coordinates, 23, 31 

Jacobi iteration, 56 
Jacobian, 25 

Kelvin, 64 

Khintchine’s constant, 69 
Kilogram, 64 
Kinematics, 31 
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Kronecker 5 , 11, 21 

Lagrange interpolation poly., 50 
Lagrange multipliers, 24 
Laplace transform, 41, 44, 47 
inversion formula, 44 
standard forms, 41-43 
Laplacian operator, 38 
Least squares approximation, 59 
Least squares regression, 61 
Legendre function, 23 
Legendre polynomial, 22 
Legendre’s equation, 23 
Limiting theorems, 44 
Limits, 24 
Line integrals, 27 
Linear equations, 12, 55 
Linearly dependent, 12 
Linearly independent, 12 
Logarithms, 1 

Maclaurin series, 5, 24 
Mathematical constants, 69 
Matrices, 11 

identity matrix, 11 
inverse matrix, 12 
null matrix, 12 
orthogonal matrix, 12 
singular matrix, 12 
symmetric matrix, 11 
trace of a matrix, 13 
transpose matrix, 11 
unit matrix, 11 
zero matrix, 12 
Mean, 60 
Mean value, 28 
Metre, 64 

Milne-Simpson method, 54 
Minkowski inequality, 9 
Minor, 10 

Modified Bessel function, 23 
Modulus, 8 
Mole, 64 

Moment of inertia, 28, 31 


standard table, 29, 30 
Multiple integration, 25 

Nabla operator, 37 
Natural logarithm, 2 
Newton’s second law, 31 
Newton-Raphson method, 52 
systems of equations, 52 
Numerical integration, 52 
Numerical odes, 54 

first order equations, 54 
higher order equations, 55 

Omega constant, 69 
Ordinary differential equations, 32 
first order equations, 32 
higher order equations, 34 
second order equations, 33 
Osborne’s rule, 5, 25 

Parabola, 6, 7 
Parabolic constant, 69 
Parallel axes theorem, 29 
Parseval’s theorem, 36 
Partial fractions, 2 
Particular integral, 33 
Pearson’s correlation coefficient, 61 
Periodic functions, 47 
Permutations, 3 
Physical constants, 67 
Plane polar coordinates, 31 
Plastic constant, 69 
Polar coordinates, 15 
cylindrical, 16, 26, 38 
plane, 15, 26, 31 
spherical, 16, 27, 39 
Polar equation, 7 
Polynomial, 1, 2 
Polynomial interpolation, 50 
Pooled variance, 60 
Population variance, 60 
Position vector, 31 
Power function, 2 
Power series, 24 

Predictor-corrector methods, 54 


73 



Probability function, 61 
Pythagoras’ theorem, 6 

Quadratic equation, 2 

Radius of convergence, 24 
Radius of curvature, 23 
Rank, 12 
Real part, 8 

Rectangular hyperbola, 6, 7 
Rectified sine wave, 49 
Regression, 61 
Remainder theorem, 2 
Residual variance, 61 
Root mean square (RMS), 28 
Roots of equation, 9 
Runge-Kutta 4th order method, 54 

Sample statistics, 60 
Sample variance, 60 
Sawtooth wave, 48 
Second, 64 

Second moment of area, 30 
Series, 24 

Sierpiriski’s constant, 69 
Sigma notation, 1 
Simpson’s rule, 52 
Sine rule, 7 
Solid of revolution, 28 
surface area, 28 
volume, 28 

Spearman’s rank correlation coef., 61 
Sphere, 7 

Square wave, 47, 48 
Standard deviation, 60 
Standard error of mean, 60 
Stationary points, 23 

more than one variable, 23 
one variable, 23 
Stehfest method, 44 
Stokes’s theorem, 40 
Straight lines, 7 

Successive over-relaxation ( SOR ), 57 

Surface integrals, 26 

Systeme International d’Unites, 64 


Basic and derived units, 65 
S. I. Prefixes, 65 

Taylor’s series, 24 
one variable, 24 
two variables, 25 
Taylor’s theorem, 24 
Trapezoidal rule, 52 
Trial function, 33 
Triangle, 7 
Triangle inequality, 9 
Triangle law for vectors, 14 
Triangular wave, 48 
Trigonometric identities, 3 
Trott’s constant, 69 

Undetermined coefficients, 33 
Unit vector, 15 
cartesian, 16 
cylindrical, 16 

Variables separable, 32 
Vector addition, 14 
Vector calculus, 39 
Vector identities, 38 
Vector products, 14 
Velocity vector, 31 
Volume integrals, 26 

Wallis’s constant, 69 
Wallis’s formulae, 21 

Z-transform, 45 
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